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ABSTRACT 


The eolar radiation reflected by targete on the erath'e 
aurface interacts with the atmosphere before being detected by the 
satellite sensors. In the visible wavelenght region, scattering by 
molecules and aerosols predominates on the interaction of the radiation 
with the atmosphere. In addition to attenuation, this interaction also 
contributes with diffuse radiation, thus changing the quality and 
quantity of radiation arriving at the satellite sensor. The purpose of 
this work is to correct the multispectral scanner (MSS) data of LANDSAT 
satellite for atmospheric interference, using L0VTRAN~3 program and the 
analytical solution of the radiative transfer equation. WVTRAN-^Z is 
baaed on a semi^empirical model and calculates the transmittance for a 
given atmospheric path for six model atmospheres or utilizes a local 
vertical atmospheric profile given by radiosound and two aerosol models. 
Thts method was used to correct the LANDSJ^ MSS data of Brasilia. It 
improved the contrast between different natural targets. In addition, the 
corrected images of two different dates were more similar ',han the original 
ones. This method was also applied to correct the LARDSAT MSS data of 
Ribeirao Preto. The corrected image gave a classification accuracy of 
sugar cane about 10% higher as compared to the original images. Some 
suggestions for improving this method for correcting LANDSAT MSS data 
are given. 
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CAPITULO I 


INTRODUCAO 

Ha aproxlmadanente cinco anos qua o Instituto da Ptsquj^ 
sas Cspaclals part1c1pb> do programa LANOSAT, e dols anos qua anallsa 
Inagans digital Izadas, atravas do IMAGE-100, sistana da procassaiaanto 
da dad^s, dasanvolvido pala Companhia "General Elatric", para extrair 1ji 
formagoas tanaticas das Imagens multiaspectrals. 

Os satalltes da sarle LANOSAT coptan Imagens muUlaspac 
trals da suparfTcle tarrastre a as transmitem para o centro da rastraa 
mento a racepgao (Elgura 1.1), onde os sinais sao convertidos am fotogra 
flas a am fitas magnaficas compatTvals com computador, qua depols sao 
anallsadas pelos usuarlos, principalmente am pesquisas qua anvolvam o 
lavantamanto da recursos naturals. 

Entratanto, esses dados coletados pelos sansoras do MSS 
nao correspondem as verdadelras caracterTstIcas das superfTcles, por 
qua a atmosfera Interage com a radiagao reflatida na suparfTcle terrcs^ 
tre, altarando a sua distrlbuigso espectral e espaclal, absorvendo-a , 
refletindo-a a espalhando-a. A nevoa umlda reduz o contraste antra os 
alvos de uma cena, difisuttando a interpretagao da imagami o aspalhaman 
to pelas partTcolas atmosfaricas conduz a uma complaxa radlstrlbulgao 
da radlagao a provoca multlaspalhamanto, acrescantando radlagao difusa 
ao campo de visao do sensor (Figure I.2),atanuando o sinal provanlante do 
solo. 


Mel boras informag5es da suparfTcle tarrastre, a dos para 
metros atmosfarlcos qua atanuam a radlagao, sao de interasse de varlas 
areas clentTficas, tals coroo: cartografia, geologla, oceanografla, agr^ 
culture, hidrologla, angenharla florastal , estudos urbanos a nirals a, 
particularmante, meteorologla. 

Alim dlsso, as Imagens servem de subsTdIo para planeja 
mantes soclo-aconomicos, levando a uma administragao mals raclonal dos 



Fig. I.l - 0 Sistema LANOSAT (Barbosa, 1975) 



■ radlancla solar direta 
Lq s radiancia difusa 

» radiancia difusa refletida na superfTcie 
Lq = radiancia da superfTcie 

Lj » radiancia intrTnseca 


SOL 



Fig. 1.2 - Radiancias Recebidas pelos Senso res Orbitals 
(LaRocca e Turner, 1975) 
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recursos naturals* o que e um dos ualores problemas da nossa epoca. 

Entre as apllca^oes praticas (Barbosa, 1975} podem-se 

menclonar: 

• classlflca^ao das cenas, quanto as suas caracterTstIcas e conH 
guracoes geologicas* tals como tipos, falhas* estruturas* do 
bras, etc.; 

> redugao de perdas na agricultura e nas reservas florestais, atr^ 
ves da rapida identificagao das causes; 

• planejamento da distrlbuigao da produgao agrTcola anual, atraves 
da estimativa das colheitas por areas; 

• melhor planejamento e desenvolvlmento urbano; 

• fflonitorla de fenomenos dinmnicos, como sedimentagao, mudangas 
litorineas, erosao, crescimentos de colheitas, nTveis das reser 
vas de agua etc.; 

• levantamentos de coberturas de nuvens. 

As corregdes dos efeitos atmosfericos nas iroagens exigem 
tecnicas que se fundamentam no model o de transferencia radiativa, para 
ca1cu1ar>se as irradiancias ascendentes e descendentes, a radiagao dif]j 
sa, a transmitancia atmosferica e outros parametros. 

LaRccca e Turner (1975) publicaram um trabalho de revi^ 
sao bibliografica contendo varias tecnicas para calculos da transmita]i 
cia atmosferica e radiancia atmosferica, reunindo-as, comparando-as e 
acrescentando*1hes algumas informagoes auxiliares. 

Os trabalhos de Chandrasekhar (1950) e Ambartsumian 
(1943) marcaram a origem da moderna teoria de transferincia radiativa; o 
primeiro deduziu um conjunto de equagoes nao lineares, de solugao exata, 
para determinagao do campo de radiagao para urea atmosfera plano-parale 
la, iluminada pela radiagio solar. Posteriormente, surgiram outros meto 
dos matematicos que tambem conduzem a solugoes exatas. Dentre elas de^ 
tacam-se: 


f 


a) metodo de Iteragao qua consiste en converter a equagao de trans 
ferincia nuroa equagao Integral, expressa por serle de Neumann; 

b) metodo de harmonico esferico qua utilize as propriedades de or 
togonalldade na equagao de transferencia para transforma>1a em 
problema de auto-valor; 

c) metodo de ordenadas discretes, qua consiste em montar urn conjun 
to de equagoes tomando discretes as varlaveis, angulos e caina 
das atmosfericas; 

d) metodo de momento, onde se atribui urn certo momento a radliji 
cia, em fungao do angulo zenital; 

e) metodo de dupllcagem, que toma. em camada sucessiva, dols valo 
res de profundi dade otica das camadas precedentes para evltar 
consumo excessivo de tempo; 

f) metodo de Nonte Carlo, que se basela em procedimento estatTsu 
CO, atrlbuindo probabllldades a partTcula com os processos fTsj^ 
cos envoi vidos. 

Esses metodos foram utlllzados para modi f1 car o tratameri 
to originalmente desenvolvido. Todos esses metodos apresentam desvanta 
gem quanto ao tempo computaclonal, adaptabilidade e praticabilldade. 

Por outro lado, como em multas aplicagoes nao sao necessa 
rios resultados altamente precisos, surgiram varies metodos aproximados 
que produzem resultados razoavelmente satisfatorios e de rapido desenvoj^ 
viroento computaclonal. Entre esses metodos vale a pena mencionar: 

a) 0 metodo de Schuster-Schwarzschild, no qual o campo de radiagao 
e suposto isotropico, separando-o em campo ascendente e descein 
dente, utilizando valores medios para esses campos. 

b) 0 metodo de Eddington, procediiaento em que se toma a media do 
campo de radiagao sobre o espago todo, sendo o campo suposto 
quase- isotropico. 

c) 0 metodo de Ramanova, que considera a fungao fase de espalhamen 
to altamente anisotropica, pelo fato de levar em conta angulos 
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Mil to poquenos. 

d) 0 metodo de Turner* que considera a funqao fast espalhanento 

como uma soma de funqao delta* considerando que o espalhamento* 

devido aos aerosols* e predomlnante na dlreqao de Incldencia da 
radla^ao. 

0 objetivo principal deste trabalho e averiguar como o 
efelto atmosferico de espalhamento Influl sobre as Imagens* para* dessa 
forma* obter Imagens corrigidas desse efelto. 

Utilizou-se no processo urn model o de transferencia radlatl^ 

va* considerando uma atmosfera espalhadora p1ano*para1e1a* homoginea* Iso 

tropica* onde a radlagao penetra (Figura 1.3) formando urn determlnado an 

gulo e com a normal* em regloes espectrals do vIsTvel e do Infraverme 

0 ““ 

Iho proximo. Para tal* escolheu-se o metodo desenvolvido por Turner e 
Spencer (1972) e* para os calculos da transmitancia atmosferica, profundi^ 
dade otica atmosferica, profundi dade otlca de Rayleigh e profundi dade otl^ 
ca do aerosol* utIlIzou-se o metodo desenvolvido por Selby e McClatchey 
(1975). 


Em 1961* Altshuler, utlllzando metodos empTrIcos, sugerlu 
que 0 calculo da transmitancia fosse baseado principalmente nas medidas 
de transmitancia obtidas em laboratories e complementado com a teorla de 
bandas de absorgao ou linhas de absorgao molecular (Goody* 1964). Selby e 
McClatchey (1975) decldiram codificar este esquema empTrIco de previsao* 
elaborando urn programa de computador denomlnado Low Transmiss Ion-3 
(UMTRM1-3), que calcula a transmitancia atmosferica entre dols nfvels, 
cobrindo o Intervale espectral de 0*25 ym a 28*5 ym (350 a 400000 cm~^) 
num passo de 5 cm~^, com resol ugao espectral de 20 cm*\ 0 calculo e fel^ 
to sobre uma atmosfera representada por 33 camadas de urn dos 6 model os 
etmnfericos: Tropical (15°N)* Latitude Medla-Verao (45®N)* Latitude 
dia-invemo (45°N), Sub-Artlco-Verao (60®N)* Sub-Artico-Invemo (45°N), 
ou m outro que utilize os dados de radlossondagem. 
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CAPTTULO II 

FUNDAMEHTACXO TEORICA DO MODELO 
2.1 - EQUACAO DE TRAWSFERENCIA RADIATIVA 

A complexidade matematlca, envolvida num modelo de transfe 
rind a radlativa, e a resol u$ao da equa^ao integro-dlferencial respective. 
Nani pul ou'Se esta equagao para determinar o campo de radiaqao no Interior 
de uma atmosfera plano-paralela, homogenea e com elementos espalhadores, 
segundo o procedimento desenvolvido por Turner (1972). 


A equagao fundamental da radiagao solar, penetrando na 
atmosfera terrestre com angulo zenital Oo(yo - cos 0o) e angulo azimutal 
expressa por: 


dl _ L(t,u,(|») 
dt 


.2'f ,1 


p(w.4».u'ti'') L(t,u'. 4'') dy'd(|i' 


4it 


“o Ej{t) P(wi4>i'-Uo*4'o) " (l‘“o) B(t) (II*1) 

4ir 


onde: 


T > profundi dade otica 

0 )^ = albedo de espalhamento simples 

B - fungao de Planck 

L - radiancia espectral total, para urn detenninado t 
p - fungao fase de espalhamento 
Ej = radiancia solar em t 


V* - 
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As malores dIfIculdadM durantt a resol ufao da equa^ao 
conslstem em determlnar as formas das duas funfoes. fun^ao do albedo de 
espalhamento simples e funqao fase de espalhamento. As aproximacoes ad 
fflitidas na resolu^ao da equa^ao (Turner, 1972) sao: 

a) Atmosfera p1ano>para1e1a, llmitada superiormente pelo vacuo e 
Inferlormente pela superfTcle; 

b) Atmosfera homogenea, Isto e, as proporgoes dos constituintes 
sao constantes com altitude, Impllcando que 


u(t) = Wo 

(II.2) 

P(t,u,4i,p',^') « p(u»4.m'»4') 

(II.3) 


c) Atmosfera Isotropica: significa que o efelto de espalhamento 
e Independente da diregao Iniclal da radlagao; 

d) Nao existe radlagao difusa penetrando pelo topo da atmosfera, 
proveniente do espago exterior; 

e) SuperfTcle e lambertlana, ou seja, perfeltamente difusora; 

f) Nao ocorre o efelto de absorgao na reglao de espalhamento; 

g) A fungao fase de espalhamento pode ser aproxlmada como: 

p(u*<l>»u' .4»') = 4itn6(w-p') 

+ 4ii(l-n} 6(u+u') 

onde: 

li - cos 4 ) > cosseno do angulo de espalhamento 

s cos cosseno do angulo Iniclal de espalhamento 
6 = fungao delta 

n = fragao da energia espalhada para f rente, depois do centre 
de interagao, onde pode-se aproxima-la medlante a equagao 


t 
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n 


Tj^ ♦ 0»95 


(II. 5) 


onde: 


3 profundidade otica total de Rayleigh 
s profundidade oilca total de Aerosol 


Para a reglao do espectro vIsTvel, a fungao de radiagao 
de Planck e desprezTvel; assim, reescrevendo a equagao (II.l), tem>se: 


w JL 

dt 




f2H fl 



PCu, L( x,ii'.4'') du'd^' 


too EjCt) p(u*t|>.-Uo»tl'o) (II. 6) 

4it 


A equagao Integro-diferenclal aclma pode ser resol vi da 
para determlnar a radlagao difusa em algum ponto da atmosfera, com as 
seguintes condigoes de contorno 


L(0,-y,4») = 0 


(II. 7) 


LjC^oiU.v) * [ I li' *♦') 
•'"1 

+ L5(To,-w'.^'l|dii'd^' 




(II.8) 


onde: 


12 - 


LjCfo.P.+i- 

LjCto.-w'i*') 
p' (m.^.“w' .♦') 


radiincta difusa solar no topo da atniosfera 
radlancia IntrTnseca da superfTcie 
radlancia difusa na superfTcie 
radiancia solar direta na superfTcie 
retletancia bidirecional da superfTcie 


A condigao de contorno dada pela equagao (II. 7) signifi 
ca que nao existe radiagao difusa penetrando pelo topo da atmosfera, e 
a equagao (II. 8) mostra quesedevem integrar as radiagoes difusa e dire 
ta, com a refletancia, sobre todo o hemisferio, centrado no alvo. 

Multi pi icando a equagao (1 1. 6) por dp e d(t> e integrando 
a de 0 ate 1 e de 0 ate 2ir e similarmente de -I ate 0 e 0 a 2ir, obtem 
se duas equagoes, respectivamente. 


d 

2ir 

uL(x,y,i)ldyd(> = 

2ir 

1 

dx 

0 ^ 

0 ■' 

0 ^ 

0 



'2tt 

'1 

'2ir 

^ Q - 


o 


p(u»i'»w' .♦*) L(t,y' »i>' )dw'd4>'dud4i 


t»o ^^(x) 


4ir 


2ir fl 




(II. 9) 
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0 

2tt 

d 


yL(x*u»^)d)idi^ ■ 


dt 

0 


0 


0 

-1 


C2ir 


f2ir 




JiSL 

4r 




'o '0 '*1 



r2n 


EjCt) 




(II. 10) 


Oefinindo 


U(t) 


'2ir fl 

L(x,ij,4i) dud((i 


( 11 . 11 ) 


L+(x) 


'2ir fo 

L(x,u, 4>) dud4> 

^0 -'-1 


(11.12) 


onde as setas f indlcam respectivamente as radia^oes ascendente 
descendente; as irradiancias ascendente e descendente sao, portanto 
dadas por: 




/2„ fl 

. 


o 


yL(x,ii,<j>) djid^ 


(11.13) 
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6*(t) • - I f dud4 (H.14) 

i© *-i 

ms qiials todas as quanttdadas sao positivas. 

SubstUuIndo (11.11). (11.12), CH.13) a (11.14) nas equa 
qoes (11*9) a (11.10). tani-sa 


^ EHt) . U(t) 

dT 


4ir 


'2ir 

1 

'21T 

0 ^ 

0 ' 

0 


P(ti»4*M't4‘) l-(T,M',a') dy'd^'dyd^ 

-1 


f2ir 


<i»ft 

4n 


Es(x) 


p(M»4»"yo*4o) 


(11.15) 


d_ EHt) . U(t) 
dx 


J2a_ 

4« 


'2ir 

0 

■2n 

^ 0 ^ 

•1 J 

0 ^ 


p(i>«4*u' >4'} LCt.vs',4') du'd4*4»d4 


-l 


JSjL. EjtO 

4ir 


2t 


ro 


P(w»4*”W<»»4o) 


(11.16) 


-I 
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Con»1d«rando a funcio fata da atpalliamnto. dada pala 
aquagio (II *4), a raprasantando n afialUlcaawnta paU aquaqia 



’2« fl 

p()i*4*p' *4') d|id# 

- o 'o 


(11.17) 


(1*n) 



P(v*4»m'* 4') 4iid4 


(11.18) 


tenesa: 


’2it ri 

P(w.4.m'. 4‘) L(T.u'i4‘) dy'da* » 4ifnL(T.w,4) 

.0 J-i 


+ 4ir(l-n) L(T."ytir+4) 


(11.19) 


Substituindo em (11.15) a (11.16) a desenvol vando*as 

obtuse: 


d£» (t) . (l-“on)U(T) . «o(l-n)L+(T) 
dt 


. “o(1-n)Ej(T) 


(11.20) 
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dE» (t) , .“oHEjCt) (H.21) 

dr 


Adinitindo-se a aproxlma^ao de Schuster- Schawarzschl Id « 
normallzando o campo de Irraoiancla difusa, ou seja 


Lh(t) ; 

Mo 


(11.22) 


As equafoes diferenclals tomam-se 


dEt (t) . (1-o»o^)E »(t) . "o(l-n)Ej(x) (11.23) 

dt Uo Mo 


dEUr) . yl-«^on) E^(T) . <»on-n) Ei(T) , (11.24) 

Uo Mo 


que ainda podem ser escritas como 


J_ [£*(,)♦£*(.)] . JL[1*“o(I-2n)] [Et(t)- E*(t)j . 
dt Uo 


»oC-2'^ Es(x) 


(11.25) 


( 11 . 26 ) 


i [Et(t) - EKtJ] . JL 0 -<«„) [E*(t) ♦ E.(t)] 

l»0 


- u E (ti 

o S' * 


Considerando 


Ft Et ♦ Et 

. s ■ I ■ 

2 


(11.27) 


Ft _ Et - Et 
2 


(11.28) 


tera-se 


^ = J_ Cl + «o(l-2n)] Ft-a»o(l-2n)E^(T) (11.29) 

dt yo 


dFt O*“o) £5(1) 

dt Uo 


(11.30) 


Como, no caso, esta*se considerando 0 espalhamento puro. 


ou seja, «Jq = 1. 


equagoes acima sao simpHficadas para 


as 
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(Ft 2 0-nlF* - 

dr Uo 


2_ Cl-tilFt - C'-inlEo* 


£i- . -Es(^l . V 

dx 


Resol vendo #s eque^oes diferencieis. obtein-se 


F+ = 


y E«e 

*^0 0 


r/y, 


o + 


Cl 




Das equagoes (11.27) e (11.28), tem-se 


eht) = K ^ _^ + _£i . 

2y^ 2 


E+(t) = (yo+2(l-n)xl _ii_ +_^ 

2% 2 


(IK31) 


(11.32) 


(11.33) 


(11.34) 


(11.35) 


(11.36) 


- la - 



V* 



! 

I 

I 


I , 
i ^ . 
[ V. 


Considerando a superfTcle lambertlana* as condladas de 
contorno r«diizaiB>se a 


E*(t) ■ 0 

(11.37) 

E*(t,) • p[E4(x^) ♦ ''“/“o] 

(11.38) 


onde as constantes Cj e C 2 sao totalmente determinadas para cada una das 
equa^oes* (11.35) e (11.36). Essas equagoes sao aquelas que representam, 
respectlvamente, os campos de irradianclas descendentes e ascendentes; 
considera>se que cada um desses e resultante de outros dois campos 
distintos* um campo calculado com p=0 e outro com pjK). Assim, obtem-se: 




V E 
o o 




(l-n)(vt) 


. r 2pp^(1-n)x 1 

° ° l+2(l-n)(l-p)T 


(11.39) 


M E 

Et(t) = 


Po+(l‘”)^o 


(l-n)(T^-x) 


♦ pPq H2(1-n)x 


1+2(l-n)(l-p)T, 


(11.40) 
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Conblnando as Irradiancias descendentes e ascendentes 
para p « 0 e p / 0» teiiHse a radiancia cm uroa profundidade otica t, 
aproximada pela seguinte expressao: 


IJ L 


+ E‘+(r)«(u+u )«(♦-♦ )] + 

® ° J 2ir 


(11.41) 


onde as irradiancias con un apostrofe representaro ocampo de radiagao pa 
ra p = 0 e 0 apostrofo dupio representa o campo com p / 0. 


Introduzindo (II. 4) e (11.41) em (II. 6), e com as condj^ 
goes do contorno (II. 7) e (II. 8), resolve-se a equagao, determinando as 
radiagoes ascendente e descendente, respectivamente 


Lpt(t»u»4) = 


4ir[y +(l-n)x 1 1 [ 


(l-n)T [p(u.(|>,y ,ir+4> ) 


+ p(y.<^.-u^,4-J] + yoP(p.4.-u^.<^o) + ^^2(l-n)(l'p)T 


-(TO-x)/yl 

.1-e J + 


r 

j(l-n) [p(y»4>.yQ.'^+4>Q) + P(p»4'»-Mjj.'t'Q)] 


8(l-n)wQ p 

U2(l-n)(l-p)t^ 


[(V^)e 


-(t+m)] 




(11.42) 


1 
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4»|ii^+(l-n)Tol 


(1-n)T^[p(-M.4.u^,.ir+4o) 




U2(l-n)(l-p)t, 


* 

(1-e - (l-n)[p(*v»4'»yQ«if+<l>Q) P(“yt4»“yo»4o)J 


8(l-n)wjp 

U2(l-n)(l-p)T, 


(ye + T-y) 


(11.43) 


Para o "topo da atmosfera", t = 0, logo 




4ir|y^+(l-n)x^ 


(l-n)T^[p(y,«|',it+<tio) 


+ p(y»4>.-yot4o)] + yoP(y.^.yo*4>o) + 


l+2(l-p)(l-n)x 


(1-e 


[p(y,(fr,y^,it+(fr^) + p(y»^.-UQ.*o)] 


8(l-ri)y^P 

l+2(l-n)(l-p)x< 




(11.44) 
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que dxprtssa a radlaqao difusa qua atinge os sensores do satellte, consj^ 
derado local Izado no "topo da atmosfera*' . 

2.2 • TRANSMITMCIA ATMQSFERICA 


A transnltincla da radiagao na atmosfera e complexa. devi^ 
do a dependincia dos coeficlentes de espal hamento e de absorgao das dife 
rentes propHedades fTsicas da atmosfera. For outro lado, a alteragao so 
frida por uma radlagao monocroroat1ca*para1ela e dada pel a lei de Beer 
(Siegel e Howe1>1971), ou seja* a atenuagao dl^ na intensidade num 
meio homogeneo, e expressa pel a seguinte equagao: 

dl = Ij^dt (11.45) 


onde: 


» intensidade da radiagao 
Y - coeficiente voluinetrico de atenuagao 
dt - comprimento da trajetoria 

Integrando-se a equagao acirna, tem-se: 

'x * 'xo t'l 1 

^ 0 


ou 




ydn J 


(11.46) 


(11.47) 


Consequentemente , obtem-se a definigao classica da tran^ 
mitancia monocromati ca 


.‘^CEDING PACii^^ 
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T . exp [- (‘ Ydi] (11.48) 

Entretanto. a equagao aclina resulta nuna Inpossibilldade 
pratica de nedl^oes de transmitancla total para uma determlnada friiquen 
da de resolu^ao Infinite* a nao ser nutn laboratorio onde a fonte Ja e 
nonocronatica. Porem, a medida da transmitancla media tav(v) num 
intervalo expectral Av e possTvel , ou seja 


T. (v) 

Av' ' 


1 

Av 


V + 


Av 


T(v)dv 

V ” 

2 


(11.49) 


onde: 


V = frequencia central 
Av * intervalo de frequencia 

No intervalo espectral cons iderado neste trabalho* vIsTvel 
e infravermelho proximo* os coeficientes de espalharoento molecular e 
coeficiente de espalhamento do aerosol sao fracaroente dependentes da 
frequencia, e a transmitancla media, neste caso, pode ser obtida pela 
interpola^ao dos valores monocromaticos (McClatchey et al* 1972), obe 
decendo a lei exponencial. 
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Para determlnar a transmltancla media total numa certa fre 
quencia devem>semu1t1p11car os valores das transmitanclas Individuals de 
cada componente atmosferlco, ou seja 

T^^(v)(total) (H 2 O) t^^(v)( 03 ) 

TA^v)(Rayle1gh) f^^(v) (aerosol) (11.51) 


A forma funclonal da transmltancla media do LOUTRAN-3 so 
bre 0 Intervalo considerado e dada m termos dos parametros Cv e u*. ou 
seja. 


T « f(C^,u.*) 


(11.52) 


onde 

Cv s coeficiente de absorgao em fungao do comprimento de onda (fre 
quencia) 

n 

( 11 . 53 ) 

u* s quantidade equivalente do absorvente 
u - concentragao do absorvente 
P(z)spressao atmosferica em z 
Pq = pressao atmosferica padrao ao nTvel do mar 
T(z)= temperatura em z 
Tq » tonperatura padrao 

n s coeficiente determinado empiricamente dos dados experimental s e 
dos valores cal cul ados no "Airforce Geophysics Lab".(Tabela II. 1). 


0)* = 0) ^ 

_P,UI 

^0 ' 



T(z) 
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TABELA II.l 

COEFICIENTES (n) DETERMINADOS EMPIRICAMENTENO "AIRFOPCE GEOPHYSICS LAB“. 


Ozonlo n * 0,4 

Vapor de agua n < 0,9 

Gases Uniformemente Misturados n « 1,75 


A atmosfera terrestre e composta de gases majoritarios e 
ininopitarlos (Tabela II. 2), quase-permanentes e variaveis, que tim partj^ 
cipagoes importantes na transferincia radiativa. A radiagao e espalhada e 
absorvida pelos gases e partTculas, pelos seguintes mecanisnos: 


a) absorgao da radiagao ultravioleta era consequencia das trans^ 
goes eletronicas de raoleculas, oxigenio atoraico, nitrogineo e 
ozonlo; 

b) absorgao da radiagao infravermelha devido a vibragao e rotagao 
de H2O, CO2 e raoleculas de O3; 

c) espalharaento da radiagao visTvel pelos aerosois e raoleculas. 

0 vapor de agua e 0 constituinte que apresenta concentre 
gao raais variavel da atmosfera; uma das maneiras de expressa-la e: 


u* = 


a(2) 


P(Z) 


0.9 


dz 


(11.54) 


Similarraente, ozonio e urn outro constituinte bastante va 
riavel, dependendo da epoca do ano e local izagao geografica e e calculado 
por: 



I 


CONPOSIC^ NORMAL DA ATMOSFERA 


CONSTITUINTE 

% POR VOLUME 


Nitroginlo 

78,09 


0x1 gen io 

20,95 


Argon io 

0,93 


0x1 do Nitroso 

0,5 X 10-** 


Neon 

1,8 X 10-3 


Hello 

5,2 X 10-** 


Kriptonlo 

1.1 X 10’** 


Xenon Io 

8.7 X 10’6 


Hidrogenlo 

0.5 X 10-** 


Netano 

2,0 X 10-3 


Ozonlo 

0 a 7x10-6 

★ 

0x1 do Sulfuroso 

0 a 0,1x10-3 

it 

Dioxido de Nitrogenlo 

0 a 0,2x10'** 

it 

Amonla 

0 a trago 

it 

Monoxldo de Carbono 

0 a trago 

it 

Agua 

0 a 7,0 

it 

Dioxido de Carbono 

0,01 a 0.1 

it 


(Segundo Farrow, 1975) 

* Constitulntes Variaveis da Atinosfera 
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w’ 


«( 2 ) 




0,4 


dz 


(11.55) 


Pode-se aflrmar que a malor parte da energla eletromagne 
tica que chega ate a superfTcle terrestre e proveniente do sol, e essa 
energla e responsavel pelos fenomenos fTsIcos que ocorrem na atmosfe 
ra. De acordo com Orunmond (197Q) mals de 99 por cento da energla solar 
esta contida no Intervalo espectral de 0,2 um a 4,0 pni. 

0 termo v da equagao (11.48) e denominado coeficiente de 
atenuagao. A sua dlmensao e o inverse do coroprimento, e expressa-se como 
a soma do coeficiente de espalhamento (k) e do coeficiente de absorgao 
(a), pelo fato dos dois processes serem independentes, portanto 

Y = a + k (11.56) 


E ainda os processes de atenuagao de Rayleigh e dos ae 
rosois sendo,tambem, independentes, escreve~se: 

® “a (11.57) 

k = k„ + k, (11.58) 

onde os indices m e £ indicam respectivamente os processes de Rayleigh 
e do aerosol. 


0 espcihamento por partTculas com tamanhos muito menores 
do que o comprimento de onda e denominado espalhamento Rayleigh (roolecu 
lar), e foi introduzido por Lord Rayleigh que tencionava explicar o azul 
do ceu. Atualmente, sabe>se que este processo e causado pela flutuagao 
da densidade na atroosfera. 
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0 coeficlente de espalbainento molecular depende da 
densidade numerica das noleculas, e a sua dependencia com o comprimento 
de or.da e aproximadaiwnte de x"**. No Lowtran-3 o calculo base1a-se nos 
resultados obtidos por Penndorf (1961) e e dado pel a seguinte equagao 

k . 9,87 X 10"^° V** (11.59) 
in 

enquanto que o coeficlente de absorgao molecular (a^) e uma fungao dos 
gases absorventes da atmosfera, da tempera tura e pressaoi a sua dependen 
cia com o comprimento de onda e altamente varlavel. No presente traba 
Iho, contudo, essa dependencia e desprezada, devido ao Interval o espe£ 
tral conslderado. 

Outro espalhamento que se considera e aquele devido as 
partTculas com tamanhos malores do que o comprimento de onda incldente. 
Estas partTculas sao denomlnadas aerosols. Geralmente define-se como sen 
do solldos ou ITquIdos em suspensao num meic gasoso, e esta definigao irt 
clui nevoa, nevoeiro, neblina, fumaga, partTculas de po» nuvens, polen e 
outros. 


Quanto aos tamanhos destas podem-se classificar, segundo 
Junge (1963), em tres categorlas. 

a) Nucleos de Aitken - lO’^ a 10"^ cm 

b) PartTculas grandes - 10*^ a 10“** cm 

c) PartTculas gigantes - malores do que lO"** cm 

Os coeficientes de espalhamento do aerosol (k^) e de ab 
sorgao do aerosol (a ) dependem exclusivamente da densidade total do ae 
rosol , da distrlbulgac de tamanhos dos aeros51s e dos Tndices de refr£ 
gao destes. 
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L(Mtrafi-3 utlllza dols model os de aerosols* digital 1zaoc> 
no proprlo prograna* correspondendo a ccndigao de ci'u limpo e unido* res 
pectlvamente 23 kn e 5km de alcance visual ao nTvel do mar* baseando-se 
nas seguintes suposigoes: 

a) A distrlbulgao do tamanho das partTculas e semelhante ao modelo 
utlllzado pelo Oelrmendjlan (1959)* Figura II. 1* ou seja 

n(r) » ar^ exp (-br)^ (11.60) 


onde 


n(r) « fungao distrlbulgao (distrlbulgao gama modificada) 

r > ralo do aerosol 

a,b*d» « constantes reals positives 

c constante Intel ra positiva 

b) A distrlbulgao de tamanhos dos aerosols permanece constante com 
altitude. 

c) As variagoes das densidades das partTculas (label a 1 1. 3) com al^ 
titude sao aquelas ajustadas pelo Me Clatchey et al (1972) cor) 
forme o modelo de atenuagao desenvolvido por Eltennan (1964). 

Outro parametro Importante no estudo da transferincia r£ 
diativa e a quantidade adimen<inna1 que mede o processo de atenuagao da 
atmosfera* denominada profundidade otica. que se define exatamente como 
0 fator exponencial da equagao (11.50). ou seJa 


T 


f Ydz 
Jh 


(11.61) 
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TABELA 1 1. 3 


MOOELO OE OISTRIBUICAO VERTICAL DO AEROSOL PARA ATMOSFERAS 


"CLARA" E "UMIDA" (McCLATCHEY et al . 1972' 
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logo para uma atmosfera Intel ra, tem-se: 

- I Ydz (11.62) 

Un valor grande de profundi dade otica significa uma grande 
atenuagao da radlagao pela atmosfera. 

Comparando as expressoes (il.50) e (11.61), pode>se rees 
crever a transmitancia atmosferica como 

= exp !“ [ ydt 1 * exp !- r sece | = 
exp 1- 7^/cose | = exp |- | (11.63) 


Na figure I I. 2 notam-se as transraitancias para diferentes 
modelos atroosfericos, calculados com o LOWTRAN-3 (Selby e McClatchey, 
1975). 


Em resumo, para se utilizar integralmente o modelo de 
Turner, e necessario conhecer a medida de turbidez, na determinagao pos 
terior da profundidade otica do aerosol e da transmitancia atmosferica. 

Entretanto, a medida de turbidez, que ainda nao foi obtida 
sobre o continente brasileiro, e urn elemento necessario, porem ausente no 
modelo. Entao, para computar os parametros necessarios utilizou'se o pro 
grama LOWTRAN, que fornece os valores de transmitancia total da atmosfera, 
transmitancia molecular e a transmitancia do aerosol, e, alem disso, este 
programa apresenta a vantagem de fornecer os parametros mencionados em 
uma condigao mais realTsta, isto e, utilizam-se osperfis verticals de 
pressao, temperature e de umidade, obtidas a parti r da radiossondagem l£ 
cal. 


RnNSi^lTTRNCE TRRNStll 
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SUSAaCTiC \V!.VT£f< 
.VlOwATITwDE WINTER 
IS62 US STANDARD — - 
SUBARCTIC SUMMER 
MIOLATITUDE SUMMER — 
TROPICAL — — 


WfiVELENGTH (i1 1 CRONS) 


Fig. II. 2 - Transmitincia Atmosferica para um Percurso Vertical do 
NTvel do Mar para o Espago a Di versos Modelos Atmosfe 
ricos. (Selby e McClatchey, 1975). 
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CAPTTULO 111 

EIIMINACAO DO EFEITO ATMQSFERICO 

0 sistema MSS (Multlspectral Scanner System) do LANOSAT 
capta os dados, focalizando alvos sobre a superfTcie terrestre em varias 
bandas espectrais que cobrem de 0.5 yo a 1.1 um (Tabela III.1). Este sU 
tema varre uma regiao de 185 km x 185 kn transversal men te a trajetoria do 
satelite, com os seis detetores de fibres oticas para cada banda. 

TABELA 1IL1 
CANAIS MSS 00 LANOSAT 


CANAL MSS 

INTERVALO 4im) 

4 

0,5 a 0,6 

5 

0,6 a 0,7 

6 

0,7 a 0,8 

7 

0,8 a 1,1 


As radiagoes totals inci dentes em cada uma das fibres otl^ 
cas sao transformadas, digitalizadas e forma tadas em uma serie contTnua 
de dados, denominados nTveis de cinza, dos campos visuais instantaneos 
que correspondem a uma area de 79 metros x 56 metros, chamados "pixel" ou 
elemento de imagem. Tais dados estao em fitas magneticas compatTveis com 
computador. 

3.1 - ALGORITMO DE CORRECT X 

A radiancia total que atinge os sensores numa profundi dade 
otica T, com angulo nadir e e azimutal 4>, pode ser equacionada como: 

Lj (t,u,^) * Ljj^ + Lqj (t,u,4i) (in.l) 


onde 
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Lj > radlancia espectral total 
Lyp s radlancia transmitida 
Lqj s radlancia difusa 

A radlancia transmitida i aquela retietida sobre o al^ 
VO, multipllcada pela transmitancia atmosferica total, Isto e 

® Lj T (T,T^,ji) (III. 2) 

C 

Entao, na equacao (Ill.l) para o caso em consldera^ao, 
deve-se tomar a profundidade otica nula, assumindo-se como topo da a^ 
mosfera o local onde se encontra o satellte, ou seja: 

L-p (o,u,<|') = Lj (r^,u,<t>) T + Ljjy (o,u,<|)) (III. 3) 

Assim sendo, para ellmlnar os efeltos atmosfericos das 
Imagens, a radlancia intrinseca de cada pixel sera determinada usando- 
-se 0 seguinte algoritino, 

Lj (t^,Wj 4’) “ (o»u,<l>)” (T q»u)j (lll'A) 


E ainda sabe-se que a radiancia intrinseca depende da 
irradiancia direta e difusa.Admitindo-se que a superfTciee lambertiana, 
pode-se aproximar (Turner et. al., 1972) como 



-T+ 




(III. 6) 


onde 


p, = reflectancia difusa efetiva do alvo 

O 

Ej^ (t^) = irradiancia difusa total descendente 
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A Equa^ao (11.44) pode ser descompactada em varias fun 
{das para fad 11 tar o seu calculo, portanto 

^DT _ 

» Fl + p 2 6 (t ,p)(III. 6 ) 


0 0 0 


0 o 


onde as fun^oes Fi,F2 e G sao respectivamente 

- — ^ " •')(t - y) + 

® ° 4u[u^ + (l-n)T^J 1 0 L 0 

♦q) + P (y.l'i-MQi'I'Q) j + y^P (y.1i.-y^j.<l>Q) + (1 -n) y Jp(u,(|.,y^,Tr + 4.^) 

T /y ^ 

+ p {y, 4 ),-y^, 4 .^) j . e y^ P (y,^.-yj,.<l'Q) | (IH.7) 


P ’2 “ 


Ztt [_u^+(l-n)T^j 


1 I 


1 + 4 u (1 -n) 


- [l + 4 (1 - n) (t^ + y)] e'V^j 


(III. 8) 


G (P. t ) = z- 

° 1 +2 (1 -n) (1 -p) T 


(III. 9) 


e, tambem para a irradiancia difusa total descendente 


■T+ 


. H [ 1 *I(tJ G (p.,J] 


onde 
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H(u tJ 2 (III. 11) 

I(t.) • 2(1-Mt. 

As fungoes fases de espalhamento da expressao sao calcuU 
das utilizando o metodo aplicado por Deirmendjian (1969) (Apendice A). 

3.2 - MCTODO DE CORRECAO DE DADOS 

Experimentou-se a corre^ao nas Imagens de Brasilia (21/06 
/74) e de Ribeirao Preto (01/07/77) procedendo os passos descritos abaixo 
e ilustrados na Figura III.l. 

Primeiro ; Nas fitas magneticas processadas em Cachoeira Paulista, 
cujas dimensoec sao de 3218 pixels na horizontal e 2352 pixels 
na vertical, delimitou-se a imagem num formato de 512 por 512 p^ 
xels. devido a liroitagao de video do IMAGE-100. 0 objetivo princi^ 
pal deste procedimento e facilitar a leitura dos pixels consU 
tuintes das imagens pelo programa de corregao atmosferica, conside 
rando que o novo formato sao quatro matrizes de 512 por 512 e gerari 
do desta maneira uma nova fita com as corregoes radiometrica e geo 
metrica incluTdas para todas as bandas do ERTS (LANOSAT). 

Segundo : Dm elemento essencial no calculo do efeito da atenuagao a^ 
mosferica, obviamente, e a representagao da caracterTstica fisica 
da propria atmosfera. A transmitancia total da atmosfera, no compr2 
men to de onda de interesse e importante e e considerado urn parame 
tro basico no presente modelo. Com LOWTRAN-3 calculou-se a transmj^ 
tancia total da atmosfera, onde verifica-se que ela e uma fungao de 
varios parametros (pressao, teraperatura, umidade relativa, visibil^ 
dade horizontal) que descrevem a condigao do estado atmosferico do 
percurso da radiancia refletida na superfTcie, assim como do compr^ 
mento de onda. Os estados atmosfericos para os quais as transmitan 
cias foram calculadas sao apresentadas nas Tabelas I I 1.2 e III. 3. 
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TABELA I I I. 2 

ESTADO ATMOSFERICO DE BRASILIA EH 21-6-1974 (RADIOSSONDAGEM) 
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Obs: Devi do ao metodo utilizado para calcular a densidade de 
vapor de agua no LOUTRAN-3, na ausincia de infonnagaoda 
umidade relative na radiossondagem, considerou-se esta 
igual a um. 
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Obs: Devido ao metodo utlllzado para calcular a denaldade de 
vapor do agya no LOMTMM-3, na ausincia do Infomacaoda 
umldado rolativa no radlossondagon* con$1dorou-so esta 
Igual a urn. 
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Obs: Devido ao mitodo utlllzado para calcular a densidade de vapor 
de agua no L0UTRAN-3t na ausencia de informagao da umidade re 
lativa na radiossondagem, considerou-se esta Igual a urn. 
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A dependincia espectral da transnltancla para con a vlsl^ 
bllldade esta plotada na Figura I II. 2. Para as diferentes bandas 
ver1f1ca-se a tendencia de aunnnto da transnltancla do comprimento 
de onda manor para malor, Isto porque a atenuagao atmosferlca e 
mals proemlnente na reglao do comprimento de onda menor do que na 
reglao de Infravermelho; e ao mesmo fato e, ainda, observado para 
a visibllldade menor. E a dependincia espectral da transnltancla em 
fungao do angulo de escaneamento esta tragada na Figura III.3» onde 
se verlfica que, para a varlagao do angulo 0^ ate 10°, os valores 
computados de transmitanclas sao quase Inalterivels. Portanto, no 
presente modelo, considerou-se o angulo de escaneamento Igual a ze 
ro. 


A Figural I I I. 4 mostra a transnltancla em fungao da pro 
fundi dade otica, onde a ultima expressa a quantidade de espalhamen 
to que ocorre na atmosfera. A dependincia espectral da profundidade 
otica de Rayleigh e da profundidade otica do aerosol estao tragadas 
na Figura III. 5. Comparando>as, verifica-se que o espalhnnento devi^ 
do aos aerosols e mals significativo no intervalo espectral conside 
rado. 

Terceiro : A conversao dos nTvels de cinza dos pixels para as radian 
cias totals recebidas pelo satelite, fol realizada pela sub-rotina 
implementada no programa de corregao atmosferlca (Programa ATA), de 
terminando-se, para cada banda espectral do LANOSAT, a radiancia mi 
xiroa e minima, usando a expressao (III. 3), ou seja, 

4i ' ‘■Ii ‘ * *-DTi (III. 13) 

onde os sub-indices 1 e 2 sao, respectivamente, calculados para 
Pa * 1 e Pa = 0* respectivamente. Assim, para cada nivel de cinza, 
tem-se a seguinte relagao, dada por 












I II. 3 - Represerta?ao das Curvas de Transmitanclas Espectrais em Fun^ao do AnguTo de Escaneaaento 







POA, POR 


1.0 




I i t I I 1 1 1 1 

* 055 0.65 0.75 0.85 0-95 

Fig. III. 5 - Dependencia Espectral da Profundidade Otica de Rayleigh e 

da Profundidade Otica de Aerosol , obtidas com os result^ 
dos do LOWTRAN-3. 


- 47 - 




255 



(III. 15) 


E a Tabela III. 4 mostra os coeficlentes angulares resuUajft 
tes das relagdes entre ci nTvels de cinzas versus as radianclas to 
tals das quatro bandas do satelite, para a imagem de Ribeirao Preto. 

TABELA III. 4 

COEFICIENTES ANGULARES DOS NtVEIS DE CINZA 
VERSUS AS RAOIANCIAS TOTAIS 


CANAL MSS 

COEFICIENTE ANGULAR 

4 

0,034344 

5 

0,042292 

6 

0,064376 

7 

0,067636 


Quatro : Sobre a fita gerada, o programa ATA efetua a leltura dos nT 
veis de cinza e calcula os novos valores correspondentes , ellminaji 
do as atenua^oes atmosfericas. Neste processo o programa efetua a 
leitura de 1.048.576 pixels, e o tempo que leva para calcular os va 
lores dos nTveis de cinzas corrigidos no 6-6700 e aproximadamente 
de 15 minutos. 

As variagoes dos nTveis de cinza computadas em comparagao 
com os nTveis anteriores sao plotadas em fungao do angulo azimutal 
na Figura I I I. 6. 

Os valores dos parametros cal cul ados na corregao da imagem 
de Ribeirao Preto sao mostrados na Tabela I I I. 5. Observa-se que a 
radiagao difusa e maior na regiao de menor comprimento de onda do 
que na maior, podendo atingir ate 55% da radiagao total, dependendo 
da condigao atmosferlca. 




Fig. Ill, 6 - Representagao da Variagao dos Ni'veis de Cinza em Diferentes 
Azi mutes. 

NCC = NTvel de Cinza Corrigido 
NCO * NTvel de Cinza Original 
ASA s Azimute do Satelite 
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TABELA III.S 

PARftCTROS UTILIZADOS NA CORRECAO DA IMAfiEM ORBITA/PONTO (178/27) 
DE RIBEIRAO PRETO (01/07/1977) 


CANAL MSS 

P 

T 

fn 

''a 

Lp^(nJW cm"^ Sr“^ um"^) 

4 

0,1145 

0,2675 

0,7990 

0,6418 

0,7007 

5 

0,1301 

0,3465 

0,8928 

0,6791 

0,5053 

6 

0,3373 

0,5640 

0,9385 

0,7123 

0,3456 

7 

0,3981 

0,5997 

0,9726 

0,7596 

0,3270 


^ s albedo medio 

T = transmi tinci a midi a 

Tj^ = transmi tanci a molecular 

Tj = transmi tanci a de aerosois 

Lpj = radiancia espectral difusa 

Nesta tabela os resultados espectrais foram cal cul ados a 
parti r das refletancias individuais dos alvos, sendo que estes aj^ 
VOS foram reconhecidos no processamento preliminar da imagem (Tabe 
la III. 6). 
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TABEU 1II.6 

REFLETUNCIAS INDIVIDUAIS DOS ALVOS RECONHECIDOS NO 
PROCESSAMENTO AUTOM«TICO PRELIMINAR DA IMAGEM 


CANAL MSS 

AGUA^ 

CANA2 

SOLO 3 

PASTO** 

MATA5 

4 

0,11 

0,15 

0,12 

0,08 


5 

0,15 

0,10 

0,27 

0,10 

0,09 

6 

0,14 

0,48 

0,34 

0,23 


7 

0,16 

0,46 

0,37 

0,37 

0,41 


1 - medida por Bentancurt (1978) 

2 - calculada por Suits e Safir (1972) 

3 “ calculada por Condi t (1970) 

4 - tornado do Meteorological Monographs, pg. 4 (1965) 

5 - tornado do Meteorological Monographs, pg. 4 (1965) 











CAPTTULO IV 


RESULTAOOS EXPERIMENTAIS 


A flm de verificar a valldade do mitodo apresentado, os 
resultados conseguidos sao apresentados nas duas formas apreclaveis* ou se 
ja, qualitativa e quantitativa, para identificar se houve altergao no is 
pecto visual e no aspecto de classificagao de alvo. 

Como primeiro passo, tomou-se a radiossondagem de Bras^ 
lia (21/06/74) e, codificando os dados de perfil vertical atmosferico, ob 
tiveram-se com LOWTRAN-3 os parametros necessaries (transmitancia, profu£ 
didade otica do aerosol, profundidade otica de Rayleigh) e, junto com a da 
ta e hora da passagem do satelite, latitude e longitude do local, acionou- 
se 0 programa ATA para corrigir a imagem original , ja antes “formatada" coji 
venientemente. 


Apos a corregao, como resultado qualitative, tirarara-se 
“slides" do vTdeo do IMAGE-100 das bandas do LANDSAT, os quais, convert^ 
dos em fotografias convene! onais, sao mostrados nas Figuras IV. 1, IV. 2 e 
IV. 3. Os mesmos correspondem as fotografias comparativas entre as imagens 
originais e as corrigidas dos canais MSS 4 e 7 do LANDSAT e a composigao 
color! da dos canais MSS 4,5 e 7. 

Nas fotografias que correspondem as imagens originais, os 
contrastes entre as configuragoes das superficies adjacentes sao reduzidas, 
dificultando a incerpretagao usual da imagem. 

Confirmou-se que, no interval o espectral utilizado, o 
efeito da atenuagao atmosferica e muito maior no canal MSS de comprimentos 
de ondas menores (0,5 ym a 0,6 ym) do que para o canal MSS de comprimentos 
de onda maiores (0,8 ym a 1,1 ym). 

As imagens corrigidas sao sensivelmente melhores; as co 
res foram realgadas e sao mais contrastantes, sendo possTvel uma melhor 
discriniinagao entre alvos diferentes, facilitando o trabalho do foto-intej^ 
prete. 



Fig. IV. 2 - Imagens Original e Corrinida do Canal MSS 7 do LANDSA 
orbita/ponto (192/23), BrasTlia (21/06/74). 


Fig. IV. 1 - Imagens Original e Corrigida do Canal MSS 4 do LANDSAT, 
orbita/ponto (192/23), BrasTlia (21/06/74). 







I Fiq. IV. 3 - Imaqens Oriqinal e Corriqida, Composigao Colorida dos Canals 

j MSS 4, 5 e 7 do LANDSAT, orbita/ponto (192/23), Brasilia 

I (21/06/74). 

i 


I" 



'3 
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Note'Se que certas areas, como estradas e vegetaqoes, apre 
sentam malor contras te am texturas nas Imagens corrigidas, aparecendo co 
mo configuragoes pouco definidas nas Imagens originals. 

Alein disso, com a corregao, a roani pul agio de Imagens de ya 
Has datas toma-se vlavel, visto que os comportamentos espectrals dos a^ 
VOS serao mals homogeneos de uma data para a outra mantendo a malor Iden 
tidade entre suas caracterTstIcas, em decorrencia da redugio dos efeltos 
atimsfiricos. Consequentemente, a corregao possiblllta a transferencia de 
comportanento espectral de uma Imagem para outra com malor seguranga do 
que as Imagens originals, nas quals os efeltos atmosfericos estao presen 
tes e acarretam comportamentos espectral s.diferentes para datas distintas. 

Na Imagem de Brasilia (LANOSAT, orbita/ponto (192/23), 21/ 
06/74) nio houve possibllldade de obter resultado quantitative devido a 
ausencia de apolo terrestre simultaneo ao sobrevoo. 

No segundo objetivo, efetuando todos os passes necessaries 
para a corregio da Imagem, tentou-se obter uma nova Imagem da Baia de Gu^ 
nabara, objetivando corregao da Imagem sobre o mar; entretanto, nao se 
obteve resultado positive porque, na Imagem Iniclal, verlficou-se que 
via nuven:; sobre a area, as quals saturavam a imagem corrigida, principal^ 
mente sobre a area continental e sobre as posigoes onde havla presenga de 
nuvens, dificultando a Interpretagao. 

Entao, para obter resultado quantitative, utilizou-se a 
regiao canavleira de Ribelrao Preto (01/07/1977), orbita/ponto (178/27). 
Na primeira tentativa, dos quatro canals MSS corrigidos do LANOSAT, foram 
ellmlnados os dels primelros canals (4 e 5), em decorrencia dos valores 
das radlagoes difusas computadas nestas faixas da regiao do espectro v1s£ 
vel. 


Neste case, as radlagoes difusas foram superestimadas, tor 
nando-se os malores contrlbuintes na radiagao total perceblda pelos senso 
res MSS do satellte, e que consequentemente mascaram outras mals impor 
tantes, isto e, as radlagies refletidas pelos alvos Imageados. 
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Isto talvez se deva ao fato de qua $e utilizou a radios 
sondagem de Sao Paulo para slmular o perfll atmosferico da reglao de R1^ 
belrao Preto, na ausencia da radlossondagem local. Obtiveraro-se com esse 
procedimento valores de transmitancia multo baixos, subestimados, acarre 
tando tambem baixos valores de radlagoes transmit Idas. 

Tals valores incoerentes devem-se provavelmente ao fato de 
que Sao Paulo e urn grande centro Industrial, tendo portanto malor concert 
tragao de aerosols de poluentes e particulados, e ainda malor quantidade 
de umidade devido ao horarlo de langamento de radlossondagem (8 horase25 
mlnutos). Haja vista que, no diagrama adlabatico, notou-se a presenga de 
duas Inversoes termicas, uma a nTvel de superfTcle, devido a resfriamento 
radlativo noturno e outra no nTvel de 700 mlllbares, devido ao sistema de 
alta pressao de grande escala domlnante. 

Entao, para soluclonar o problema, procurou-se uma radios^ 
sondagem continental mals proxima possTvel de R1 belrao Preto na mesma d£ 
ta de passagem. So fol possTvel a de BrasTlia. 

Repetiram-se todos os passos e obtiveram-se valores de r^ 
diagdes difusas menores do que os anterlores; mesmo assim, nao se conse 
guiu realgar o canal MSS4. Isto porque a Informagao dIsponTvel sobre a ca 
na de acugar neste canal de Imagem original e escassa. (Tabela IV.l). Pa 
ra todas as fases de crescimento de cana de acugar consideradas, no inter^ 
valo de nTvels de cinza entre o superior e o Inferior sao encontrados ape 
nas 2 ou 3 nTvels de cinza. Portanto, Imagina-se que ausencia do canal MSS 
nao afetara significativamente a classificagao desse alvo da cena. 

Finalmente, para comparar ambas as imagens, original e cor 
rigidas de uma mesma data, quanto ao aspecto quantitativo, tomou-se como 
referincia uma area delimitada, totalmente coberta por fotografias aereas, 
para servir de area teste para medir a precisao de avaliagao da area cana 
Vieira, obtida com a Interpretagao automatica da imagem original do LANOSAT 
e da Imagem corn g1 da pelo programa ATA. 


r..«li:irtrM DAS DIFERENTES FASES HE CULTURA W CAHA CfUMnO 




I 

I 
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N« fotografla airaa obtave-se ima classlflca^ao da cana da 
acucar aUborada palos patquisadoras do Projato ”EstatT$t1cas AgrTcolas" 
(EAGRI)* do Dapartananto da Sansorlamanto Ramoto (DSR), do Instituto da 
Pasquisas Etpaclals (INPE). E. para o computo dasta araa, ut111zou-$a uma 
malha da pontot aquIdUtantas da 1*0 mn* obtando* asslm* uma araa corras 
pondanta da 106*03 km^ da cana da a^ucar numa araa total dallmitada da 
197*21 kiR^ da fotografla aaraa. 

Os rasultados comparativos da avalla^ao da cana da a^ucar 
antra as Imagans a as fotograflas aaraas sao mostrados na Tabala IV. 2. 

Em ralaqio as fotograflas aaraas* qua sao consideradas co 
mo s^do vardada absoluta* obtava*sa com a Imagam original do LANOSAT a 
classificagio da 68X da araa corraspondante a cana enquanto qua utilizan 
do-se a imagam corriglda obteva*sa 77%* tambem am ralagao as fotograflas 
aaraas. Devido a dificuldada oparacional do IMAGE-100* nao sa conseguiu 
delimitar a mesma araa total nas imagans do LANOSAT* mas obtava-sa uma 
araa total* quasa similar* da 180*33 km^. Isto corrasponda a urn arro da 
8*5% na datarroinagao da araa total das imagans do LANOSAT utilizando o 
atual sistema oparacional do IMAGE-100. Obviamanta* as percantagans da 
classificagoas da cana mencionadas seriam maioras para ambas as imagans se 
as araas totals* isto a* a da fotografia aaraa a a do LANOSAT, fossem 
iguais. 

Anal isando os “slides" das imagans obtidas por equipamento 
DICOMED (sistema fotografico do IMAGE-100 qua transfere os sinais corra^ 
pondentes dos pixels diretamente da memoria para o material fotografico), 
verificou-se qua axiste uma distingao radical na classificagao da cana de 
agucarentre a imagam original a a corriglda. No primairo, notou-se qua 
sao classificadas aquelas canas qua sao fortemente iluminadas. Logo, para 
a imagam original axiste uma tendencia de classificar alvos mais ilumina 
dos do qua na corriglda* a nesta foram classificadas aquelas canas menos 
iluminadas. 
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TABEU IV. 2 

l^REAS COMPARATIVAS DA AVALIACXO DE CANA DE ACOCAR ENTRE AS FOTfl 
6RAFIAS ACBEAS. A IWGEM ORIGINAL t A IMAGEH CORRIGIDADO UWDSAT 


MfkTERIAL DE 
CLASSIFICACXO 

XREA TOTAL 
DELIMITAOA 
(km2) 

AREA CANAVIEIRA 
RECONHECIOA 

(km2) 

OlFERENCA 

REUTIVA 

(*) 

Fotografla 

acrea 

197.21 

106.03 


Imagen 

original 

180.33 

i 

72.23 1 68.1 

1 

Imagen 
corri g1 da 

180.33 

81.98 

i 

77,3 
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Para avallar a$ precisoes da classlflca^des da ambas as 
Imagans* amp11ou*sao*s11da"quecont«n as duas Imagans, slnultanaaiaanta, 
com urn projetor KODAK ate atingir o tamanho da foto aaraa e fez-$a a com 
paragao. 


Nasta comparagao* obsarvou-sa qua cartas araas qua nao cor 
respondlam a cana foram classificadas come araas canavlalras. Isto aconte 
ceu davido ao fato qua certos alvos tern refletancia aspectral samel hante 
a da cana, sando pertanto, uma deficlincia do sistama am si. Csta ocorran 
da inplica radugao da area real correspondante a cana; porem, como as 
ta fato ocorre am ambas as imagans, original a corrigida, acradita-se qua 
os rasultados relatives, isto a, a comparagio antra as imagans original a 
corrigida, nao sejam alterados significativainente. 

Portanto, para a imagmn corrigida, no computo final da i 
reas, teve>se urn acrescimo da cerca da 10% da area canavieira classifies 
da. 


Porem, uma conclusao nuantitativa mais segura, expressa em 
termos numericos, somente pode ser obtida com urn experimento bam mais ala 
borado a controlado, a partir dos princTpios estabelecidos para este meto 
do de corregao atmosferica para imagans MSS dos satelites da serie LANOSAT. 


CAPlTULO V 


CONCLUSOES 


Demonstrou'se » aqui> que o presente roetodo de apllcagao da 
equagao de transferencia radlativa nas corregoes dos efeitos atmosfericos 
nas Imagens, produz resultados relevantes. 

Pel os resultados aU *'gados, observa-se um aspecto visual 
inelhor e um aumento quantitativo no reconhecimento dos alvos; no caso t^ 
pico estudadOt a regiao canavieira de Ribeirao Preto, obteve*se um acre^ 
cimo aproximado de 10% na area corretamente classificada como cana de agu 
car. 


Os calculos tim sido coerentes, embora nas faixas de com 
primento de onda menor (0,5 •* 0,6 um e 0,6 - 0,7 um), a concordancia seja 
menor (Sharma, 1976). Isto em decorrencia do espalhamento atmosferico ser 
roaior nessas faixas, como tambem a resposta dos alvos, principalmente a 
vegetagao, ser menor nessa faixa do visTvel do que nas faixas do infra- 
vermelho proximo. As analises mostram que as Incertezas no conhecimento 
exato dos valores dos parametros amblentais, utilizados como entrada no 
presente metodo, podem resultar em erro:> ^istematicos nos calculos. Por 
tanto, devem-se tomar precaugoes nas especificagoes dos parametros sobre 
os quais o metodo esta baseado. 


Alguns aspectos podem ser implementados, tais como: 

a) Absorgao pel os aerosols, usualmente dispensavel na regiao espe£ 
tral visTvel, mas na regiao do infravermelho proximo pode ter 
alguma importancia; 

b) aproximagao da superfTcie como lambertiana, comum, por^ existem 
muitos casos em que a suposigao e falha; 


c) formulagac bi ou tridimensional para o metodo, na corregao da 
imagem devido a presanga das nuvens que alteram a radiagao pr£ 

veniente dos alvos; 


-rjj 
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t 

I 

j d) descrl^oes mals exatas das condlfoes atmosfericas locals, com d^ 

I dos medidos slmuUaneamente com a obten^ao das Imagens; 

i 

I e) utlllzagao de urn programs classificador automatico SISTEMA MAXVER 

I (classificador Gaussiano por Maxima Similaridade) mais sofistic^ 

i do do que SINGLE-CELL, no IMAGE-100. 

} 

I 

I A inclusao dos refinamentos acima mencionados podera pro 

porcionar urn metodo de corre^ao atmosferica para transferencia radiativa 
mais sofisticado e urn si sterna operacional de corregao atmosferica com con 
digoes adequadas, produzindo melhores resultados em fungao da realidade 
terrestre. 

Por^, nunca deixara de existir o problema intrfnseco do 
si sterna: a semelhanga espectral de alvos, que sempre limitara os resulta 
dos alcangados, por mais sofisticados que sejam os model os de transferee* 
cia radiativa. 
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DETERMINA.CAO DA FUNCAO FASE DE ESPALHAMENTO 

Considerando a teorla de espalhainento da radia^ao eletro 
magnetica pel as esferas dieletricas homogeneas, isto e« espalhamento de 
Nie, obtein-se o$ parametros necessaries. 

A expresao analTt 1 ca» que e fungao do comprimento deonda, 
parimetro do tamanho e Tndice de refragao das partTculas, define-se por 


P(COS x) = ^ . (|Si|2 + |S2|2) 

2nx2Q^(m,x) 


A. I 


onde 


= fator de eficiencia de espalhamento 
m = Tndice de refragao 

X = parimetro do tamanho 

r = raio da partTcula 

X = comprimento de onda 

Si *82 = amplitudes do espalhamento 

Os termos acima definidos sao dados pel as seguintes reU 

goes: 

ur^Q (m.x) = r (m.x) = r ( 2 t+l)(la |2 + [b J 2 ) 

^ ^ 2ff t=1 X, z 


A.2 
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10 • IHj • il02 


A.3 


Sj ■ r a n, ♦ b,t, 

1=1 1(1 + 1 ) ^ ^ ^ 


A.4 


$ 2=1 — ■ |a,T, + b,n, I 

1=1 1(1 + 1 ) ^ ^ ^ ^ 


A.S 


X = 


2irr 


A.6 


Nos trabalhos de Stratton (1941), Von de Hulst (1957) e 
Deirmenfjian (1969) pode-se encontrar urn tratamento mais rigoroso. 

Na equagao (A.2) a^^^ e b^ sao coefi dentes de Nie dados por 


I i 

_ Vg (mx) Vc (x) - my,, (mx) %(x) 
c(x) - m\ 


, • 

_ mv, (tnx) (x) - 'to (rox) y, (x) 

io'i'^{nu() <;(x) - \{mx) ^l(x) 


onde e q sao as fungoes Riccati-Bessel , e os apostrofos indicam dife 
renciais era relagao ao argumento da fungao. As fungoes e sao dadas 
por 



1 1 


“ 

r 


ff 
% ► 



I 

1 
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n^(cos x) 


dPj(cos x) 
d(cos x) 


A. 9 


T^(COS x) 


cos xn^(cos x) - sen x 


d|nj^(cos x)l 
d cos x 


A.IO 


onde P|^(cos x) e o pollnomlo de Legendre. 

Determlnaram-se os pararoetros aciina e calculou-se a fun 
gao fase de espalhamento para as partTculas de aerosols de ralos compree^ 
didos entre 0»1 e 5»0 pm com Tndice de refragao media de m^ « 1,33. 


A Figure A.1 inostra o resultado calculado, onde se verifi^ 
ca a distribuigao angular da radiagao espalhada pelos aerosois, para dife 
rentes comprimentos de ondas. Pode>se notar por comparagao com a Figura 
A. 2, que o espalhamento dos aerosois e totalmente distinto do espa1hame£ 
to Rayleigh. 
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Fig. A. 2 > Oependencia Angular da Fungao Fase de Cspal hamento 
(La Rocca, 1975). 
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16 IF n.ri.Nl) f.w(o, I 

IF (Hi .(,r.Z( t) ) JPT 

* IF (K’M PRIM a3«, I » ? . I )• (CH(K. I ). K = 1 . 10 ) 

fH<V»r)=(H(O,|)*l,0 

17 C0*.Tl"t 

170 IF (IrTl.n.CO. 1 ) OP TO 9 

IPs-; 

I K s f* 

XleHl 

CAil ' CIf.T (Hi . T» , r PI . TA. IP ) 

JU' 

T X 1 S. T V ( d 

ro :r- ► s j , <■ 

It) r ( *^ ' = ' V ( I- ) 

IF ( ) ’•vt f .t • . n r-r r'> 

1 F M ’ ' - r . F ^ . .1 ) H s 7 ( ' I ) 

IF (f M F .f.i .9. .( 1 St* T.j 
19 if (r‘.Lr.M.^ .f. . r>.- rM .oT.u) jl*jl^l 
J?rl I 

IF ( ] T ''►‘f ,f . . j 5 r.f. TP 

rci I '•'M ( I f . f . 1 <. ir ) 

vVr* 

ir c* • ^ 0 * ) j f • ^ 

?u Cl : . 

r h f - . t f ( X 

IF » • • -'-r 


71 
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CO*'Ti‘ If 

ir (jt.fc.j?) t*i»txi*vn-£mi«» ji ) 
c*«** Niw ’'rfUE VC«TICaL ‘'AIM (lUAhTlTirS Vh(l-tt) 
ir (jo.tc.O) ►niAT <*?o 
00 

>l«Z( T ) 

X2»7( T*n 
If 

If (I.fJ.J?) »?«»-? 

If (l.f-.Nl.) u 2 » 7 ( r 
OS . 07 

(••••• vfh'‘'» T'^ajECTORY 

PXr{Kr«»l)/(ot*X 2 ) 

THE Tt-tSl N( SPhI )/C A 
PHIstSi .( bPHT..iX)/C4 
f)ET*T»iMA-PHt ! 

sat 

If (S^^'l .I.T .1 .£- 10 ) D 5 s(Rf** 2 )*SlN(bn*Ci)/SPHl 
fETa*5r la^^rT. 

PS?=P£Ta«P'*I-ASiLE 

pH|.l jn.-pwi 
SRsSR *?3 

00 ?a -< = 1.'' 

EYsOS*' T ) 

If (l.'-.\L) r.J TO ’? 

" If (£-!(•<* I ) .f'J.J , 0 . )R.frt(A,Ul ). EC. C ,0 ) CO To ?3 
If (E I <. 1*1 ) ) 60 TO 2 « 

EYsOs*( T )'£►*(■<• I ♦!)) / Al-CSf t'-C'* I !♦!) ) 

.GC TC ■ 

72 IF I ).r j. J.O) GO TO 23 

• If (E-f 1*1 >.£w.e..)T 50 Tu 23 
If (?-)(•<» I ..rj, £*-{<. I-l j ) 60 TC 2 a 
FY=f.V/ 'i.O i(f u(K, I-l )/ru(K» I ) ) 

GO TO 5-. 

23 * EVsJ. 

2 « VH( K)sV-( ^ )-r V 

If (J’.CS.O) P-INT 433, I ,X i.( Vf'(L ) .L« i.ti) »PSI .PHI ,HETa .ImETa ,3* 
IF (1. . jL) 6T TO 25 

If (I.1.E>J9) £H(»,I»n.TN 
If (I.r..Jl) ?H(4,t).r»l 

RNsfHf J, I .1 )/t- (R. I ) 

SPhIs-.^hI ."X/- . 

If (SlLP. jf-.o-*) S5»HIr>3wP 

25 CO'.'TIvjf. 

GO TO 1/ . • 

C»»** Hjol7''.fAw PAT- 

26 DJ 2^ •'*! . T 

t»(Os >4 '•£*.■«{•<,)) 

If (P I'll . *.T, j ) *(^ ), (• .Gl* rx(-. ) 

?T CO’TI lir. 

GC TP «■« 

2« CO.M •: 

* 1 **^ -••J 

K 2 s «3 

If .f 1 ) 

J2r M . ' 

TS*-Y 
J».l’ ♦ 


I 

I 
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If (N3t.e4«l.*Nj.H{.'ilv.I( GO TO JU 
CALL 30iNT l»r2>rK*-«*l«>*7*TA*IP> 

DO 29 

a* W(0«r»(K) 

TX?»T<(0 
VN7»V I 

If (H»,LT,H») HmH? 

J2»i< 

IF TX?»T<l*v„2-Cn(9.N) 

IF (H».;T.hi> T<1.T<2 . 

IF ( J! .CO. J?. A ,a.H?.LT.rtl ) fNl»T*2 

30 Ao»(<<**^i 

IF 

CO 31 T«l.Jl 

I )-BE 

IF l-H£ 

JUI-JeJ 

IF (HM ^3 ro 32 

31 C0MI;lJt 

32 XaM.M . ! ■ 

IF (H a '(.Lf..«) io n 3 a 

call JIlLf (Y'Y.ON. jP.Ta.1P) 

TX1=T<f9) 

IF(TX l.LT.O.O) TX3«TX(9) 

•IF(TXJ.lT.O.o) TX3«TXC>) 

IF ( Ji .i0a;.4.\fD. r.ui ) GJ rc 33 
HMT'.s '.'»/! A 3-9E 

IF (A 4S(x-h-,t.% ), r.T.'3.3 '' j 1 ) an to 3 ? 

33 IF ( jl.to. *•'.»%:. va»TX3 

IF ( J1 ,UF, J2) tA2«T*3 

IF (H’.'.C. Ht) Ta?«T<3 
'I' ) J2«,a 

IF (H’.'JC.' .n-i.l-2.LT,w aIn J H»N«lI., 
IF(H?.lT.h» »v) j?s j 
IF (Hi.LT,» I.) PRi^iT a<J0. HMI\ 

CO TO 30 
3« , Pilr.T «3X, 

I f (WJ.CT.-l I j TO 33 

IF (T TYPE. ; 3. JP.H?, jr.Hl ) PrtlNT A3F 

ITypf ,? 

Tx? = E jf ), I ) 

Jrtf-.= I 
J2*l 
P2» ). ) 

N * 0 . P 

c*»«* N"' ''■FKf vC-)TICiL PATH CUAnTITIFS Yh(l-fc) 

35 IF ( J’.CO.O) P-ip T J?)- 

00 40 

P£f sCu( j) 

IF p.fstni 

IF •4t.^•Tl\P 

IF *“,'') •«»F«TY< 

If ( 1 . • . t ) Y 1 5 i ( J ♦ 1 ) 

X2s2( i> 

• IF (J." •. J'' .* «T j A^»M 

IF (J.t . J-I‘ . - r . x:.s'“ . I) »?8|-.'IN 

IF {• 1 . M'. , r. I , A, ) A?*hf^ 
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THrT«^SlN(S®M| )/C» 

S4LP«.X»S»*H1 

If (A»s(*j-HM).iT.i.0E-b) *tP«*sn,rsAtP}/c* 

. PET»AtP*T«£ T • 

If ($3>JI.GT.1 .OE'IO) 0<a(nC*X2)*Spi(Ptr*CA)/SPMt 

THCTAal jO,C-Tm£T 

BETA* -.rTA*?? T 

PS I »!* ; T A" 41 AnoLE ♦ 1 0» 

SR«S^ *■*5 
00 30 

AJs^. h(K»J) 

If (J.fl.jl) 

If (J.f *j»«(K) 

If <j,rc. 

If ( J, • * . J^‘ I »•. Ar.r , 4RS( ) ,LT . 1 ,CE“5 ) «J*TX(K) 

If (KJ.rc.C) 00 TO Tf 
If (j.f’J.j?) bj*^(<) : 

If (J.f-. Aj*T<(<) 

36 IF (f j.i C.C.n.p^.aj.rQ.y.w) f,y ip 36 

If (ij.EP.FjT OJ Til 3f 
FV*f S.(Mj-fcj)/4Lr f,( AJ/Pj) 

CO TO M . 

3f EV»r£.’.j 

• GC TO •>? 

38 FVr(.' 

3^ VH(r.)r' H(K)*f V 

If (Jo.tr.C) P-iI>T «35. J.X1.(VH(L),L»1.6).PSI.4 LP.(JET*.THETa,SH 
If ( J.^j. J^.tf.C.R?.f,r.Wl ) ijJ TC 65 
If (J.r.. j‘'I»’.a»;0.k?.F 0.1) (it) TO <i? 

If CJ.’E.n PlvaWf F/rM( 0 ,J-l ) 

If (J.T.,J?»1) litn^rr/TA? 

If (J f .. J?. ANC.HZ.ro.O «f*PLf/Y»,P 

If (J.‘’..(J>'lN*l ).^M0,K^.tU•l 1 MR*®Ff/TA3 
If ( ? • Lf . fit , P\ ) Px«i,o 
,SPv-Ts':/Lf.W*i 

If t J. r .. JP . t\0. k P.fS.T'i tiU TO “1 
40 fCKll ' E 

Al IF (»<• »'..l£ .0 ) lO Tn «T 
If U .( (,.( T ;>rtIM <*38 
IF (I .fO. C.T '5'J r-j u7 
IF (I P' .fr). 1 ) HWJ*. I a?" 

K 2 *I 

Tl**? 

If (/• ).L, .n.nji ) uO TO t»r 

H*k-l I 

IF ( f (J • 1 1 J a J- 1 

R r ^ r T 

PHc Ip . SPHl )/CA 

TSs''^! 

PSs* 5 ’ 

u c r • * 1 ^ 

0 ? F ( K U i I ► ) 

r.u ’f 
n3 PF^ - 

PS’-; , " 1-* ^ 

5 h r r . . - I • 

C - If* 
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PMI»Pm 

DO 0* • 

4« VHCK)p?.'‘VH(ir)-E(K) 

.60 TP «7 

«s 00 

«6 VHrK)s2.C«VH(K) 

PCT*«?.J«bET<» 

SP»2.r«SP 

IF (Mj.Ed.Hl > 00 TO or 

nhmi%\/yhi 

sPMi*<;p.( 4dr,i r*CA) 

IF (SbmI.lT.BN) SpHTaSPHl/MS 
■ 60 TO I*} • 

•F COsTItUt 

PR|t.T 4J6. HM ■ • • 

DO AP K«l.lo 

«6 C0*T1-Iit 

«9 BSiiTF (6»*19> 

PRUr 0*a2l) (;mI ! 

I«1 
L» J. 

I VI «V 1/5,0 

IV7aV5/5.*.90 
lVl»5.TVt . 

IV?»?>»TV2 

IF ( lwl.Lf.3^1.) TVl«T5ft 
IF ( 1 v7*r.T*5pOO'J ) Itf?«5y0o0 
— ■ IF (ny.LT.S.) DV»5. 

•*" I0var)>> . 

lV»lVl-IPv 

lCfU»T=U 

C****' BFOI* I -f. Or THAKSMiTTANCt calculations 

50 IVsIV*TUV 

• lF(JP.»t.o) r,o 10 52 

IF ( I COUdT.t 0, G ) (,fi TO 51 
IF (ICTJ^ I*C0.5u) 'I. I TH 51 
CO Tfi S2 

51 ’ ‘ lCnUM = G 

PRIM VJ 

52 PC 5 3 

TX('<)a'i.O 

IF (K.l. (.«) TXOOsl.l* 

53 CO'JTl.'.'L 
ICni!NT= ICOUM^l 
SU*i»o," 

v«iv 

I«( I v-35r )/«>*i 

IF ( r>.i T. i'»oo) r,n TO PI 
IF ( l".LT.?7/m) r,o TO Po 

"'.LcCUlA^ SC’.TTf.RjNj 

CBt-o,-'. rr“?o*i / 

SUf'sS ;"*Ta( pi 

IF (! .•.lT.qjol) .r.ft TO -c 

IF (I'.LT.I.TdoO) CT TO o9 

C/U-4r 

IF (Iv.'.*'.J?3«00) r,3 TO 3« 

IF ( ; J.CF.^rsu J) r,(j To 5S 
00 TO hT 
5A XXaZOr.'J 
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xi»( V-1JIM'J,0)/AX*1 ,0 

Ll«l 

L2-S3 

GO n Se 

5S XXaSn.i.u • 

XI«(V.^7SoO,rt)/»X»Sr,0 

Lla'^r 

LP-IO-* 

9* 00 sr 

X0sXt-rLn4T<:<> 

ir (X n sv!i'».9r 

5^ CO'jTl '!£ 

58 TX(:»). U^)«C«(N) 

CO ▼I iO 

58 TXfi)«i:-*(rt)*xO*(CH( O-CSCN^m 
• TXf ^)s in)*TX( 1) 

60 SU*'«S' '*Tx(<*> 

ir { IV, -,f, 14500) io n 47 

GO Til <.4 

(•••*«**• »4T*:4 V-P0U4 CO’JlIiXUJH 

61 ir lO n 7? 

IF (U'.tT.ToO) jO TO *>•> 

X|«( V-7‘!(». )/SO.*l, 

no M 'iP»i*)s 

XHmXI-'’L04T(»iH) 
rr (X-) 

63 COfiTT -ME 

64 TX(‘*)r'*bf :,S) 
on TP 

TXf ')»'•%{. '<H)*XH»(CS(;W>-C>r<H-l ) ) 

60 TO »i7 

66 TX(5)»f -(“oTo. )• 

67 T*(5)-i(S)*TX(3) 

SO>i*Sii»<*T X ( ® » 

00 TC ’2 

^••****** '<iTH.“,"ni COnTP'U''’’ 

66 IF (I.,.LT.?i ‘0) GO TO ^2 ’ 

^K4sI*1/«6 

5 J»«»5 »'«♦ T X ( 4 
fO TP 7 ^ 

. iirw \ *PDu^^ 

69 IF ( I . .LT . ) 4'iJ. I V. -jf. Vrt75) GO TO 70 

IF (1» .if .l4'5?O.A<,J,Tv.ljr:.li'*uO) GO TO /\ 

GO TP '6 

70 

GO TP * . 

71 

72 Kl.) 

IF ( V f ) ). L T. 1 , ;E-7 . , -n To 76 
*‘SI =Al'\,U.( , f 1 ) )*v; U T ) 

IF (•• « 1 .11 . - Ut A ) n I.) 76 
!F (T . ( 1 . T. A? ) ,n ▼o /5 
IF ( . I I . : , .1 ) ► 1 a 

ro TT . rl ) ,<•» 

IF <r 5.1 .1 ► .f ►.(<)) Gl TP 7 a 
t3 CL> '1 • . 

7« TX<11-T. («)»nA(K-,,-jr^ll..j).(fn(K)--St)/(F«(K)-FW(A-l)) 

or 7- 

Tb TXf 1 )t . 

7t- Cf ’ I - I 


I 


♦ (. f'lif.o r,»ScS 

IF n-.LT.>C*t.AKir.T v.r.t.bOu) CO TO 7T 
IF ( Iv.LT. Ul«0*«NU. Itf *bT. U9Fo) GO TO Fb 
CO TO fll 

CO TO 

j«(Tv-ir‘>5n)/5*i?i'i 
IFf .i(?)*LT*1.0f?0) r,n TO «J 
Kl.l 

^))*CF( J) 

IF r,n 10 8J 

IF (t<s2.r,|0,I»6rt?) Tj o< 

IF (>*«?. fit. 2,0) fcl»*0 
00 On «K 1*67 
IF (»•«?. LC*rM(<\>) 00 T'' Si 

0 CU».TI i:f. 

1 TXf^)3TK(^)♦(T‘i(K-^-T«>lK))•(M•(K)-»s^)/(^*»lK )-fnck-i ) > 
CO TC *'3 

2 TX(2;*0.0 

3 CO»;TI.ijF. 

• ••*•••* i7:it . • 

IF ( I •.LT.^7«i,C«.I v.OT. J2FCI CC TO AT 

L«T-6s 

XU! 

IF ( V-'( ^ . t ^ • 1 • Jt”?« ) f,0 TO «T 

IF (vi^’.L I .677« ) »io nt t,r 
IF (t-;0.r.t.;4.'n4^) r,0 •»u on 
If(tii'».0r.l,*>)Klaj6 

00 •*« ' sK 1 *nT 

IF (*s'**l-k.'Ff’»'v)) <i 1 T'^ 

4 CO»TI I'E . 

5 TX( T ),Tt<(iv )'f T«iK )-Tf’{''“l 1 )*(FC(K )-wS3 )/(FOlK )*FO(K“l ) ) 
GO TP OF 

16 TXfOls'.d 

iT COnTi i;L 

;•••**««• A,£irSi;i C aT 1 *. c T i "!i. 

4LA*'si .t,r *c/v ■ ' . 

XXsf'.’. 

IF .ni.l.) f.J Tf! VO 

00 I' ^ *< *! * *^ •; 

XD**-L ( '. > 

IF( X'' An 

«8 CCfTl-i.t 

89 yxs(f7f: )*C7r'-.-i))« <r>/f v*u. )-v>(s-ii)*c7(*<) 

VYr(( *f (*. )-' 7AC.-1 ) rfA(M-*XC -l ))»C7*l»!) 

90 ix( If. <7 ) 

TX{7)-!<i*a(7> 

so* *f- ' *TX t 7) 

TXfTl.A-' 

0 0 '*‘••10 

IF (Y TP v<; 

IF (T .f*. ).UF.</.1 ) .» T-' '»1 
IF (T.r^ ..) .1 T-^l ti 

1 A{ •»'! - 1 V M I ) 

f.G •i" 

CC 'P '*'• 

Y Xi • 


2 
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«1 T)l(>O«0. 

f« COnTIMIC • 

T«(in)»1.0*TV(l0) 

TX(9).?X( 1 )*Ta(^)»TV( 3)«TX(9) 

IF ( I .« .tii: .130U0 ) T<(3)arx(a) 

IFf jP.r.i..<>T*(9}«T*(T) 

ABal.'TXfV) 

IFf lV,«-.i.IVl,09.IV.ro.f «^2) AM.C.5**0 
SU*<AbS'i^<a*ap»0'/ 

|F( JP.c-O.o) *'K(Tr('>,a2?) 1 V.AL*m»T»{9).(Tx(R)»'<« 1»F).T*( io)»SUMA 
IF (Iv.CtC.IVP) UO TO 9S 
60 TO so 

95 P£*C IXY 

AB«l.r“i'JHX/CV2“V1 ) 

PRIM u2fl. IVl . IV2. SU-A* *0 
PHiM -lor.lYy 
IFdXv.tC.O) GC Tf, too 
GO TO (96»?.nT.96,l00)*lXI 
9* REaH -OO. V1,v2.Pv 

AV*t»l '0^.0./ VI ' 

ALA'l>* /JCV»/V2 

• PRIST olft VI • v2»0V> ALArit Av<< 

SUMA«so 
GO TP 49 

9T IF{nrr‘’L.EQ.n) oO T’J 200 

GO TU 

^8 REaO '4 fiO»MOPEL.lHA2F,lTTPt.LEI..JP.Tl»,Mj,H2.H3,HL.H0 

PRT'.T -u-o# '’ontL* Tha I T Y^-I^LEs* JP» IM» Ml » m 2» ‘i3»ML* 90 
~i GO TO AOO 

too STmf 

. Aoo FO-«'A'(iniJ.no.3) 

• 01 rC»'*AT ('’Etc, 3) 

• 02 F0»“AT (T*). t. 3»FA, i.pLly. 1) ) 

• 03 . FOO 'At C«lfc.3*2rT.« ) ) 

•o« ro4"A» fisr^.p) 

•05 ros’'Af toto,?) . 

•06 format (Trio. 3) 

• 07 .rC9‘*AT (//10»«?lM .4^017001X1 . PAIh» ALTMudE «*r7,3»llM nM»RASGt ®* 

1F7,3.1“ <<M) 

• 08 rOe;*.tT (//1Jv.50h SlAnT PaT.» UFluEPf; -LTITuoES hi Ai*0 «;» *Mt«C Hi 

I»»r7,).5M h' «»PY,l.lflp ■lM>7tNlT4. Ai« ,Li »*F7, 3 *Hh DEtiKEf.S) 

• 09 FO-'wAf (//TCv.JaH S'.A*|T P<T^ to spACE Mo** AlTITUDf ‘<1 *.f7.3.l9H 

1KM» Jf'JlTh A»'ftL€ JEu4t£S) 

• 10 rtj4«AT ('?■><. HH MjOrL mTHOSPhE'T Ml'llH * T‘'C*’1 Ca'-) 

• ll pOJ>‘4t (/?>9ii,Hh k'j.'CL AThCoPhF'mE tII>2lH s MlOU A 1 I f o'JT SU“hEP ) 

• 12 fO“‘'At (/24x,1-h Ali,CSPnEME •Il*flH * MlULATIfCDt hINTtk) 

«1 3 rCa'AT (/i'0>,l^n i'jTL AT«Ci>PhF»’F ,Il»eii « SUi-AHCflC SUhM ij ) 

• la rO>»;'AT (/?tU,l4.. wjiel «TrUiPHF*-E »ll*<iH » 1962 fS S'A’-Oa^O ) 

• 15 ro9”4T 1/20X.HH « I'PL ATM^iPhri'E .lt*?lH * ^Ijt'-AKCTIC ..IsUP ) 

• 16 FO = "At (/70«,l-»- ma/i ftri-L *Il«3r- . »a5»13h VISUAL > 

• 17 fc4''At •'nfi »".rb,i." vlsu'*s pANr.t »i sla 

• In FCs"A» (/if-A./iH r^cTt.r^CY xAf.rt vi« pf/.i.liH ch-i to v 2« •r/.lti 

lAH CH-! foR f>v *,f^. 1 ,Ch ch-1 (.Fj. 2." • ••.F5.?." ’<ICr(Jt.S )") 

• I 9 FC’i'At l/lPv.j*-- F.i'KAvFf.I StA LLVIL AfS(;«crP j-hIJo ,tS// 21 v I 1 (im<.ai 

Itr VA;*'‘1> Cc-2 f’C. L7‘'SE MThG.,EK (CO'.!) m 20 (Cf"'T) 

2 MTt SCAT .•rt.''SnL LiCPf ( 0" X ) Z/*- *• 7Hf,n C A •? » 1 0 A • Tma*- • 1 

30*.CH•.T^ L“. 1 0* • At . '«» * 7 hi,m C» “2. U X*?-r.n. 1 3X*2 Hkm. io»» ahA ’.m t») 

• 20 FCC'iAT (1>-1./// Sj.M VcRII^AL M*('FJLtS t>A T. "PS I "* 0 » • "PH I t > » •• 

lPtTA",'»*."T*riif,<i<,"i. At,r,t"J 

• 2l FOW'A’ (/iPA.tf" M l-'l)»0(L14.J)/) 

• ?2 FC"”"A’ M -1 . / I P A. .T4,< r«E^ "AylLEf.MH TIMAL m? w t 5 ««hCPz ♦ . i . < 
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, • 

l«H020*r. n 2 C0N7 CJNT HOL iC*T »C«OSOL ACROSOl IMTCftHATfO 

2 Cm-I ►*IC'*n»'5*8(*jc^-iTHANS).«*.20H *hS ABSORPTION ) 

• 23 rO«M*A» <luX*T6*lOFV.').ri?.?> 

•2* rCarAT (« iNTtcwAieo aSoRPTION FR0P**»15a" T0»»I5»* Cn* 1 ■•»rio»2» 

1- .AVfnftOr Tmi,Srt!TTANCr ■".rA.«) 

• 23 FOR*'*T (1u>t»^Flu.3) 

• 2* rCPMAt (//0)i,«»£POSr.L St*lTLRI».a NOT CCfPUTtOt |MA2F«0") 

• 2f FORMAT ( lhl.///l0»./>0M rt0hl20NT»l PPCFlLfS/) 

• 28 FliPtAT (KX." *.l««.rT.T*"jkH*M<«’'.F>'.3»"AM.«NOLE««»F8,A»''GtOM, HtUr, 

IF «-.tr,r.’*<»P.SET»»"»rA.5^". VlS*".F6.n 
•29 FCPPATf 3F10.T.2F3. |.?rlu.3*2FlC*3) 

• 30 FOBMATf »C),»«TnPuT -<F TtOHOtOalCAL 0* T A I "/ J OX • ■•»FT • 2» « KM» Pb-.FF 

J,2.« • a. TB*'.r5. 1 C» Oi." P I . IF''P"» F5 . 1 • " g» HCL HU>< 1 0* T YB*»Fb . 1 . 

2- t* w?:j lC»<ITT«**.lPt.o.?,«.(;>, ozone OEnSI 7v«"rE9,2. • c 

jTSI'aL B*t.v/F«".0pr6.l»" ^K•fiA* CE»"*F 10. 3." kN " ) 

• 3| FC>MAT(<*(f A,p,;r7.5)) 

• 32 rOn» A’ (" STAfTINr, fABAi-EFCHS AND Ai.GLE MAvE BEEN HI OE FI NE D J >il « 

1 ".FI. ,3."Af.r.LE ■".no.o) 

• 33 FOP»'A» (" TMjfCTOS<v MTiS£S eAaTmS *TmOSPMJ.Ke.- CLOSeST DISTANCE Of 

j APPP-SCH IS". F 10.2. I (. /. 1 X. "Enu OF C-LCULATIOn") 

• 3a FlircAT ( litX. t4.F6, 1 . 1 1 ( (.lu. 3) > 

• 35 fcp''at (n.rfl.i.Ari'i.3."Fv.4,f7.i) 

• 3a fo^iiat (• » ".riP.J) 

• 37 format (» PATH ir.Tct'SECTS t^RTh - PATH CHAN;.£0 TO TTPt ? nUh m? • 

1 0..1 

• 38 fOR'A’ (" Ctir>ICt OF T.i'’ PaThS FOR ThIS CASE "SHORTEST pATh TAhEN. 

\ FPN !/»f r.^R pbTm Set 

• 39 FONi'At (« choice of TbO PnTrtS FOR "mIS CASE "LONGEST PATH TAkEN. 

1 fin ^jort ctih SFT LE'i « 0 ") 

• •O F3ai.*T I" H? SrT l rbS THAN ovl*' ANU HAS HfCN HfSfT [00*0 TO 

1 Hf*I I * «l • h 2 » ".Flo. 3 ) 

• •1 Fllci'ATC* *ir,r,L AT loSPUtPE hO. '"./ A <."2 ( kh)". 3x.**P <MB)".aj(» 

j "1 (T) PiH f tPH H;-''(r,.i.H-3) c3(3M.*'"3) mo. OEN.") 

• •2 FOr»MAr(" FPr CO'.riTlOHS hAT Exist AT SEA EwEL FOR ThIS vISuAL hi 

IKCF"./.- T{ i: Tnt-1 ASsUilt The TrAvshITTanCE OUE to FOG IS r.lvf x 
28 t THr TRmNS" IAnCl *'' O.bS MjCPO'«S") 

STnp 

El.n 
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SU^^C'fT ( «> VH»M»NP»Tk • (P ) 

CO**i*r j /( J«»»P( 4 » »'I|.T(^. JP )•£►»( 10» 34 ) »!•>»«» «L» PC* Cp»Co» PI 
OI^fNSlJK TXMO) 

SUo^niiTUi pnl^T COIPOTCS ThC PCF<4>*Tlvt INOCI AflOyC «NU B^LOn 

A r,Iv"J AcnTyit *<t) I*tTEKPJLAT£s r*pu»aMl*LCt TO OEICHmInE TmE 
EiJijlVM^M AnSl^PTH AMnu*Jrs AT IhaT altitude. 


X IS T-*»’ oysTlty^ 

• Tx(9, TmP mCA.« i<rm«CT»vE I.jQICEs AdOwE and rLLOm k 

N IS T-*C LCvtL iNTrr.r*’ CJ=*HESPUKoTMi Tu i oP ThE lEvEl BElO»i i 
% p ■( If * COnCI ffS OTm hOU£L ATpOSPHt«t LCVEL .If NOT MP ■ 0 
T»(l-«» 4PP *Huu-<Ti PCP »<M *T height X 


N«VL 

NP«0 

IF (X.LT. •).■!> <«0. 

IF (x.r,r*.^( <L) ) on to i 

00 1 I « I » yU 
w*t 

IF (X.?( I )) 9 , 4 , \ 

1 CO*;^I •!" 

J?»K 

FAC«( <•/('<))/{ /(.I? )-’(►<> ) < 

PXisPf *(P( ..J?)/P( ri.n) )**f AC 
T»i^T'***V)«("( *»j?)/T{‘-'»H))**eac 
V(* l« ►.'<••» *)«f 4 >)«*FAC 

TX(?)jC . j?)/T(". j?i-4,i6L*P«i»‘^(H. j^)«I(M. j2)*-C4 

Y>»0. T < I )♦! <( J) 1 

IF '.1 TO ) 

I) J T « « i . n 
T X ( S ) •« '■ . ? 

If ,).r-, , j.o) *,C fO 3 

IF (€ O.r, 1 .K. )J.0> in TO J 

T <f Otf j^T/EH(<.n))»»FaC 

3 CO'.''! 

GO TO -J 

4 L'«I 

I*'M'>.iC.n)G';Tis 

0 : ' « 1 . ' 

5 T 4 ( •: ) t" 'I ' »•! 1 

4 T f ( n ) ,r I . 

Y • . 1 • ) 

• c*..*. ri(. ", p Pt 4.^ •) J Th»OiJuH 5*« 4Sl SO Lt>S0l4 HC.j.iMEO 

1' )-i. J 

« c f 1 1 i: 

TF •'■•p.T »00. a»N*NP. I *( ■* ). Tn. IP»( T «( n) .K*! ••» ) 

T<( »)r’ '( . 

t'i.Y* .1 . 

PtT'jF , 

OUO f'J'»"Ar l/." FHC ' I’JIMl Hl I liH Ii" , rl 0,4 t " KN, 4 «* » I 3 , " » NH«" » I <T. f 

1. »T.rY \.*Ovr \ ('fL H l’**"*.l3*F»IPX."r fUlV. 

PAM'IM* '*■. *’i” AT 

r , • 


• •'•.‘‘aU.D 
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SUano<|T!N£ ( Hi .H7* a'K.Lc * dl •L£^>ML ) 

•PI*'tH;|}N TamO) 


Th!^ ;|i|iJ*U(tT|hC rAi.C<JL*TCb ThC XnITIAl 2£nITh *mclC‘ ( AN ftCE ) 
T**IW; INTO icCJUhl oEThACT IO.>i AffCcTS alwE'< hI.h?. anO bETA 
(■Mtnr pn* li The fA«TH CC’^TSr ANOUE SChTENOEO at Ml ANO M^ )• 
ASAimri'b TMr rCEAACTIvE lu-JtX TO CO-oSTANT |s a r.IvEN lAtEa. 
TOO fiuEATEH ACCoFACt IMt»tAiC THE Nt'MdtH OF LtvlUi IN ThE mOoEL 
AT»l'5>’Mr.NE* 

This ^u^-HouTtSf ca» hcmOwed each th£ ppcoRAt* ir nct hcquIaeo, 


ASP (.)>AAS1AI A) 

ACOS(? )»AhCC^( 4) 

IP»A9 

CAsMl/lhO. 

Jtl»Pf*Kl 

k?>Ar«M? 

IT.<. 

ei.Pl.CA 

IffPl.fC.O.O) 3l«AtJS( F//*l ) 

TA*c»v?»Sif;(B»)/(A-<«ci'StH»)-xi) 

THET«tT«N(TiHG> 

ir ( T* F1 .lt. «. u) TMtT^THri.HI 

spwi .'p.nHf Tj 

Asr,«T.M/CA 

c’ PPPT •!>)«. f |.ANr,.TANr, 

TA.Tm; T 
TH»p.- .*>»CA 
J ANf.l te’» ET • 

rpi«ft. 
fETA*-. 

PtTl., 

Pfl7»- 

rpTj . 

FHT?. .C 

If (‘>1 .1.1 .«•< 1 -iJ T ) ? 
c P^t?T <ivr. IT 

Y 1 

If f»0 TO ^ 

P l P C. 1( r) jp T- A 
I" J . T. CI'S(M j-M 
ir <y p “•-1 ) c.'. ,« 

2 HPf «k 5 

H 1 -Jf** . 

3 AK ',1 r b.** J 

Iprli -cl; 

5Pt-l«’ . . 

AK.a.’ '\.f /(.t 

C pm* T i <*. »i . A »f,. JI--M 

A IPelO 

CAl L . i M ( U) . f‘ , . . <K. ! » . I •') 

Jl*‘- 

T * I » I * f ■» ) 


ORIGINAL TAGi: » 
OF POOR 


C*LL a' 3 HT ( « 2 . t K. Np •<<». 1 K» I P) 

JinK 

ir (No.EO.n stN-i 
ir (Jt.EO.J?) T<l«TXl*V>“tH(9.Jl) 

00 y 

j) 

X 2 tSE.»(J*t ) 

IF (J.Fw.JJ) 

IF (J.i4.j2) X?xBctH 2 

S*LP«<t*SpHi/<? 

«LPai<-T'p(S*L°) 

RNxEHf 9»J*1 )/E*i(9. J) 

IF ((J*l ).F 1 .J?) 0 4 «rN/£M(<J. J) 

IF (J.E3.J>) h\sEH(^.J* 1)/T*1 
IF (( J*l ).*:o. J.Tj. Jj) RKsYH/TXl 

BET*T4ET-4LP 
Far-TiN(ALP) 

IF (J.Vt.jl) f 1 sFB*TAn(TW£T ) 

F9T»F -TpFr} 

BETAsiF T4 *^ft 
THi=TyrT/CA 
BEa'lET/Cft 
C*1LP/CA 

P.RfIT <02. J,7( J).T><FT*ALP»HET*ri£T»*FaT*FB»THl*BE»C 
IF (Xa.CO.^t »H?) CxPl-AL® 

IF (SiLP.OP.PS) RN«1 . 

SPMla^iLP/'^M 
THi.-Ts'.'lN(SPM| ) 

COVTI- IE ■ . 

IF( U.LE.O*0) r,j TO 29 

60 TO ?o 

COsTKiIl 

TAf:F,s-TAf;<, 

Aft(,Lj:-s'’I“A'lOI t 

ANr.sft-.'rUF/CA 

PHT*(T aCA. *n <6.14^6 

IF (Ml.ur.O.O) r,n TO 3 

CC'.TI IJE 

IPxlOi 

call o'’If!T (Ml . rL.N.MPi »Tx* IP ) 

TXIsTvf 3) 

VM=Y<- 

IF (LPl .Eg. 1 T 

ir (»'.pr,.7) j.r*HL 
Jlsl! 

IF (!••’. '•£. *11 > .;a TO 13 

CAI L -"IM (t<2. rl.,‘j.*JP. I <» IP) 

TXPsT 9) 

Y».??Y . 

J2s». 

If (Jl.EC.J2l Ta?sT‘,J*TX( >)-rn(9.J1 ) 

XlrlF«7(J*l) 

X2rPF*7(J) 

If (J.Fj.j!) Alsr.L*^-i 
IF (J,t*...i?) A 

SAl Cr.l /K? 

1 -n ' » • 1 •'(pY'i-f'F 


n 
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c PRIM «U?. J.kl*Z(J|>S»Hl»S«LP>HMIN.RC 

ir (SaLR^lT.I . u) GO TO it 

SAlPr<;PHI 

ir (M-IN.gT.mJ) CO TO to 

11 ALP»*ST'.(S*LP) 

THfT.iSn (SPhI ) 
l’jT«*LP-‘rhCT 
PtTl««FTl♦^fT 
re»TA» f iLp) 

IF (J.‘E.jl) F0«Fa-T*Nf THtT) 

• FBTl»rcTi*r« 

THtBl-rT/cF- 

PCeFtT/CA 

AL»ALP/CA 

• C PRINT uO?, J.K’.TH€T.»Lr*»tfi;Tl»H£T.RMlh.HMlH»FHTl»THl»BE»AL 

IF (XP.FO.Rf *h2) c«»'I‘AtP 
REr*r»fV.j) ! 

1C (J.PC.jl) Kf.Fry.,t 
IF (J.FJ.J?) 

IF (J.r j.i) r.n TO i? 

PN*EH(0* j 

IF (J.C-..JJ) - ssYM / rw( 9, j-i ) 

1F( J.pP. J?*l ) Pf.rPEf /TV^ 
lFfJ,r0.j;<) t F/YNk' 

IF (F'.LP.or.t'N) Rn»I. 

~ SPM»«;»LP*Nf. 

IF GO rn u‘ 

r.c Tc iG 

12 XIrJif 

iF ( AJStZtJ'-'’2).LT.l.0t-lO.ANC.J.ME.n Go To 13 
CT TP 1“ 

13 J-J-1 

' XUar*?(J*l) • ■ 

IF (J.C'.jl) *l«.‘E*Hl 

IF _( j.n. j?. Ai.i! '.ur.jl) *l*9t*>*2 

14 Xifr'fr.7(J) 

HPIi s,l«r.M'!-»»F 
IF .r. j; G'] Tf) 25 

IF (7f GT TP Id 

Pi r»f -''-* 0 ) 

IF ( J.r J<) firFey. 

SAl P«/1 /J<2 

Atf .f S**l P ) 

Thf > 

HFi=.ir-ThCT 

THsTA f^lH)-TA.( TMC r > • • ‘ 

n«T .'srP T? 
f.ci'r I ;?♦(•! T 

pt 1 f TU'iLT< 

Al.t^l T/CA 
• TMsT„rT/i..' 

C PHI' T a ,1, *2*ThfT I r-.o*. T2,F »_T ,PMI^l ,hHIl. »FrtT ?, TH I »H£ * AL 

f N.f » f /» M '<..1- I ) 

• If ( 5>iL^ • (it . ) P'iaJ.O 

SPmT=^ tLP.Pf 

CC ■'«’ 13 

IF T X-«« 1' 1 *T; t f(9, jl ) 

Y M r 1 * 3 

IF ( A‘-5(F i-"?/ J« I ) ). LP. 1 . Of 5 ) Yf.l*TX(>) 

IF ( .• .5(1 f : r J* n ) .(.F . 1 .Of 5) Yf.i,TX(9) 

P N r * » 


60 Tfl 1<» 

CALL pf’PfT (MM|n,ym,n,*IP,TX,IP) 

7X3«Tv<9) 

IF (J.r:.Jl.x\C.H?.r,t,wi) f.O 1C IF 

IF (j,r<;.jl tfii. j.rc.j^) TX3»Y|,2»TX(9 )*tn(9. J) 

ir T' TxJ«T*m 

IF ( J.‘ 'i.jl .X* ) CL TO XT 

RN«rrr/rxi 

IF ( 5 '.'-P.r.r.eN) PM» 1 . 

SPwl»5eLF*R'i 

x■xl♦s"^I-'’r 

CIFbA ) 

IF (r;r-l.0E-5) 19.l9.lo 
X2«Ff .“VIK 

T»-rls^Mf.(SPMl ) . • 

IF(Pf 1 •i'1 Fin3B-TA».'(TriEn 

If (r .tP.!) UZ Tn V(i 

CNts(»ti-i.0'>«ALrr,(rT)(’*‘1.0)/t‘'‘tF-1.0))/(X?*Xl) 

F0T3=-Ti* ( TiifT >•( I ,0-1 ,w/l 1 ,0» TX3/f X?*u*.X ) ) ) 

PETsO.'itFl-TPtT 

BET? = :«-T2*‘‘ET 

pm'.c !"T:^orTi 

IF ^ <J j TO ?3 

PtT=H; Tl*< ,*ot r? 

09l*P*-l:n 1 

C3?=HfT’-f:i 

r tJ 3 . J A « ( f M t n 

IF(rp-..Ll .Lr 1 .A .u’Pi) '*,0 TO 25 

IF(nr.?.cT.r»,^) ;,c rr ?p 
iF(r,f '-.f.T.ruj ) r.n to ’ti 

pf T A e-4 r 

Faisf •Ti*i>.c*(roT?*FPT3) 

LE*=J. 

&C T'J -iC 
ftMr, Cl*r.pi;- 

F5TsP r;*f P»?»FrT3 

Prtl'T 'iil, J.-FTA ,Pi<T,rcU ,FR12,FpT3,Txl ,ymI 

CC TP ’t 

PET .f ♦ ( » f ^ 1 ‘if T? ) 

LC* '1. 

F0Ts?.'’*(f*’T1*PvT2»r»TY) 

PKT’T •••,.1. J.-E U.^PT.PoTl »» HT2«ffjT3. Till »Y»%1 

IF j.» f .^-i 1 g: n 

IPrJO 1 

IF f>->'.rr,.n 

sp,-:s iT ( t' .-I i ) 

IF (7l 11*1 ),| l.nr) '.0 Y;; «« 

P N r T > • / V * J 

IF (< -I.r.Y.f N) FK*1. 

SPl-lr^P'-I /•-• 

THfTr 'I*V-Pnl, 

Cn TP •' 

C»i I •'■I* I (*•;.»> . ..•p.ix,l>‘) 

1 m ri t Y»«, i.( ji ) 

F f. r ; 

J2*.l 

If 1 .l.f .1 » f • »1 . 
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THrT«\?I» (SPuI ) 

r,c Tn f 

DCTAaurn 

len«o. 

FBTaFjT i 

?6 THcTai* ;»Lt*(0l-.5FTA )/(l .♦f mT/TANG } 

O0PT4*''tT*/r* 

B«nr.Ti /CA 
THl«T^tT/C* 

PHT'»T f r T A.OPtTA.Pd T» THl. T-KG 

IF (T-rf.r,T.TN. jP.T.<rT.cT.TH) T^ft >( TN«Th )/d . 

TMjsT-fl/C^ 

PflT'T .».•*, frcT I »P.F «T,Tnl 
TM aT /C A 
TMJaT ./CA 

PR1' T .1 S, T**,T '.T-Jl .T»'l 

spu]«?r o«(T) • 

TAf ;,«T.\ .f T*^FT ) 

ITalT«.1 

nHr=Aa‘i(‘n-»4rT4). 

OTA;sA^S(ANr.LF-rhET) 

IF fKrTsO.*'»(*"<iL£4ThtT ) 

IF (I’.to.K.) nj n 

IF (f> .CI.1 .oF*F.a iO.nrH.iiT.l.uE-7) au TO I 
28 ANrLCaTnf T/CA 

— ■ PRI'.T >*'.»■, A»i(,Ut*ir 
RETI.R 
29 Hle‘<2 

ANr.LCaf/CA 

• PHJI.T •u'!. A»'^lF,IT 

RCTI.P-, 

C 

AOO FOPI’At (//" TITr.PATinN •.Urteit.R ••-*13»//) 

AU] FlJt ( I'l'f lo.F.hF j 3,' ) 

402 FCA^AT ) 

403 Fo'fAT (" » JiF s"f I 'i . 6 • " PR»".E16.‘0 

404 fUF'AT (•• Ti TAL («t fA 3 1 1 « . .'> » F 15.6 * "• FUT » •.t'ltt.o." IHFT s"»F10 

1 .6."T i* = 1 U.6) 

405 FL»’-AT (?Fl?.f.) 

406 fCnAr (3*./i «7,;inH A iI,lF 3"»F^.3." DEGHEES i RtCUMt-OTEf) 

1 FpC*. 1 .E a i„i f IT f^ .<A1 lUN". I3.-* )" ) 

in’' 
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apEndice c 

PROGRAMA ATA 


PROGRAMA ATA destlna-se a: 

a) Ler os dados locals ' 

OA - numero do dia do ano* em contagem contTnua | 

GMT - hora media Geenwich 
OG - longitude (grau) 

ON • longitude (minuto) 

AG > latitude (grau) 

AM - latitude (minuto) 

b) Ler os nTveis de cinzas originais da imagem do LANDSAT 
NIV - nTvel de cinza original 

c) Calcular as radiagoes difusas. 

d) Imprimir os nTveis de cinzas de 5 linhas, com radiagoes difusa, 
intrTnseca e total, nTvel de cinza original e nTvel de cin 
za corrigida. 

e) Gravar a imagem em fita magnetica com nTveis de cinzas corrigi^ 
das. 


[ 


riLt i»*'CPi .ut*iT»T»t>t.pf:coRn«a6 

riLE 5 • A* Vt 

riLC 6«TSU."M I»T«PE,Rf C0»0»o6 


OlH:KMC;f. C(^J2),mM?),MloW).Cir5l2).NC(5l2)»Nl»(5l?).PR(«).p*( 

DATA cR/.y*R,,«J?i»..Qj'»5,,.>(f?e/»PA/.6<4id..679l.,7i?3.. ^596/* T«A/, ? 

• 67S».J’lh^)..5AH..S9^7/.iL/.ll*ia».13'>1..3323*.39Sl/ 

CPsTm ; •», 

REtO(‘.lf’t) f>A.aPT.or,.nM,Ai.4H 
100 ron'-AK^r 5. 1 » 

c**«* CO'-vF-'i-O CJrBj, 00 Sus-SATll, Em vi«Au 
KL1 

lFfi’r,,i.T.u. ) g: TD 1 , • 

W»^L1 ’ 

r -0 TP •> 

1 h»-»L! 

FI1«Aj5('Ag)*An/oO, 

IFf ir,.LT,o. ) GC TO 3 

Firll 

CO TC ■* 

3 F«-Fn 

A f 1 f »F «n. i) 1 7 

CAl L O'^LAh (PA«G^T>wnN«f IA*^mO*A2C>2SO>asO) 

C***.*SCl FC-? f't ‘TakCA 

rc lA ti=).h 

CCTpo.6»7,*)I^ 

5 C0*?Asr..S3 

CO ’’0 ’ 

« C0f"'is:-.>5 

GO TO 'i 

7 CC'I'Aal.Tb 
GO Tp ^ 

8 CO»."As 

9 CaiL r.'.iC f7»'J./»’A.ft«;'^,Mri,PF<!.C0N''lA) 

PCT«-;l.t.C(T;>j ( : I) ) 

P 0 » =- M cc-c P'-'f 1 n ) 

POs*“ f f !•})) 

TPi^ SsT./i ( I..> 

. ALPt>ii.(!c») 

n i*( - '.'♦I,. «S,.1 ,A)/tPu>'*Pi.A> 

e‘.r*t,-PrA 

C«**«*CA1 c. ■'•*.5' Ci;rf. pr c*'*:’'. 

cat L "**.f 7’ '.r^t. ♦A’ I'f 1 »P» «;.1 ‘-ANS.mu A ) 

C*»*««Lf ? Tit ; A ; 0 'u:u P).' ri'2A 

no !T 

()»’'. '•'tj*. r f, )*i|v 

1 0 f ’.a ' ' 'y ' 1 ^ 

I'f If 

If f I 1 Cl . a>. j.i.T , !▼ I ) ir 13 

If f ’ . • 1 / ! . • i. J. if . t ’■<•<*• >.J.I-T. 3a? ) ^0 Tw 12 

If ( ' . I ’ . ’ C 1 . * .. j. , f . 3/ j > .. j ir n 

|f('» tti ’••t*,'»<'*c*** t,^,L3,|71) iiO T^ 13 

Iff*# Gil If) 12 



\l 


ORKUNAL RAGE 18 
OK POOR QUAUTV 

lf( uO TO U 

ir( !,*,t.1bI.*n 3, j.U. 1^1) vjO TC 13 
ir(I.'‘.T.3«l.»r.3.J..iT,1"0.Bi»0.j.i.T.Ta2) wO TO 1? 
irrI..*.T.3U.»(v3.J..JT,i«l) 00 TQ 99 

It 

»NrA«iLFft(2) 

PI-PFK?) 

CO TO 1 *. 

12 A9>0. 

*NFA« ’ '.P*( 1 ) 

Pl.PFiM) 

P2.rr->(1) ' 

PT«f'l.'*2 
CC TO 1“ 

13 *S«a.»l 2 <l 
*NrT«iL^*(3) 
pi.pri(3) 

P2«Pf?(J) 

PTsPUSf 

C****C»L'‘I)L>' 04$ F|iNC-rS ri . F< £ r, 

U rl*(l./( ‘ • ( 7 ’ i*£'^t*rOT ) ) )• (E‘j£« (PJT-ZMA )«?T*Z <C‘p2*ENE*2m»*p 

*'*T«TrA* 9-Z»i'’.n?.TP4 .S' 

r2.((Z •:*/ ®CT ' ))*i» . ♦“••Z‘•A•t^£•(l. ♦‘•.•FnE«(F 

• 0T*Z«a > )«T^A»'S ) 

/ I t .*? . *t'<F •( I ) *PQT ) 

PCsC 0 «'“*(ri*r 2 * 3 ) 

C( j>»- 

P(J>»4*lrA.r(J) 

PIf j)-»« P( j)-«0- '/Tt'A'iS 

CIj .•)/4*Fr. 

IFlCIc ci( J)s«.I( J)*!>.5 

j).L .c( j)»j. 
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